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ABSTRACT 

It  is  recognized  that  a  means  of  converting 
electromagnetic  fields  into  visible  images  would  greatly 
simplify  the  problem  of  designing  structures  to  optimize 
their  compatibility  with  such  fields. 

Two  needs  must  be  met  to  allow  such  visualization. 
First,  the  fields  must  be  transduced  into  a  phenomenon  which 
is  capable  of  being  imaged.  Second,  there  must  exist  a 
sensor  which  senses  that  phenomenon  and  generates  an  image 
visible  to  the  human  eye. 

The  bulk  of  this  work  was  to  investigate  the  parameters 
and  limitations  of  one  possible  imaging  procedure,  called 
Dye  Transfer  Modulation  (DTM) . 

On  the  basis  of  that  work  and  the  development  of  a  new 
sensor,  the  AGA  Thermovision  device  was  recognized  as  an 
exciting  alternative  to  DIM.  A  brief  investigation  of  its 
potential  was  conducted  and  is  presented  herein. 
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I-   SIGNIFICANCE  AND  APPLICATION 

A.  THE  PROBLEM 

The  design  and  analysis  of  devices  which  interact  with 
electromagnetic  fields  is  of  principal  interest  in  many 
technical  areas.  One  such  area  directly  applicable  to  the 
Navy  is  that  of  shipboard  antenna  design  and  placement.  A 
second  is  the  problem  of  countering  the  coupling  between  a 
structure  and  a  strong,  broad-band  signal  such  as  an 
electromagnetic  pulse. 

In  both  cases  theoretical  approaches,  including 
computer  analyses,  become  severely  complicated  by  the 
existence  of  structural  discontinuities  (junctions)  and 
assymmetric  configurations  of  the  structures  under  study. 
Even  simple  models  of  such  structures  show  significant 
inconsistencies  with  experimental  measurements.1 

It  would  be  highly  convenient  if  there  existed  some 
simple  and  direct  means  of  measuring  the  currents  and 
charges  on  such  structures.  One  candidate  for  such  a 
process  was  suggested  by  Keigo  Iizuka.2  A  more  recent 
investigation  of  process  was  made  by  James  Dodd3  to 
determine  its  application  to  the  problem  of  structures  with 
junctions. 

B.  THESIS 

This  thesis  is  to  determine  the  parameters  of  that 
process,   to   characterize   those    parameters   and   their 


lR.  W.  Burton,  The  Cross ed-Dip_ole  Structure  of  Aircraft 
ill  fL£  Electromagnetic  Pulse  Environment 

2Keigo  Iizuka,  A  Method  for  Photographing  Microwave 
wiin  i*  Polaroid  Film,  Technical  Report  no.  558,  division  of 
Engineering  ancl  Applied  Physics,  Harvard  University, 
Cambridge,  Massachusetts,  March  1968,  p.  3. 

3Dodi,  James  L.,  A  Direct  Method  of  Visualizing  the 
Surface  Distribution  on  Antennas  ajE  MTcrqwave  Frequencies 
wIITi  IppIicaTion"  1°  ^r2§sZ."^i£e  Structures,  "Piasters  Tnesis, 
U.  5".   Naval  Posfgraduale  School,  flonfereyj  Calif.,  1974. 
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interactions,  and  to  determine  a  set  of  parameters  which 
would  optimize  the  application  of  the  process  to  the  problem 
of  imaging  the  scattered  fields  of  structures  with  junctions 
and  therety  to  the  determination  of  the  currents  and  charges 
present. 

C.    THESIS  OBJECTIVE 

The  objective  of  this  work  is  to  optimize  a  method  for 
imaging  electromagnetic  fields.  This  method  is  called  the 
Dye  Transfer  Modulation  Process.  Optimization  is  considered 
to  have  occurred  when  the  images  obtained  with  the  process 
show  the  best  possible  definition  of  the  fields.  Ideally 
one  should  be  able  to  easily  distinguish  regions  of  low 
energy  density  from  those  of  high  energy  density  and 
determine  the  boundaries  between  them.  Further  the  range  of 
energy  levels  distinguishable  by  the  process  should  match 
the  range  of  energy  present  in  the  fields  being  measured. 
In  short,  it  is  considered  necessary  to  find  the  settings  of 
process  parameters  which  would  maximize  the  dynamic  ran^e  of 
the  process  and  match  it  to  that  of  the  fields  being 
measured. 
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II.   FUNDAMENTAL_BAS|S 

There  are  three  fundamental  bases  to  understanding 
dye-transfer  modulation  and  to  explaining  why  it  may  be 
successful  in  imaging  a  field  when  conventional  photography 
is  not.  These  bases  are  (1)  understanding  the  basic 
photographic  process  (2)  understanding  the  spectral 
distribution  of  energy  and  the  application  thereof  to  the 
basic  photographic  process  and  (3)  understanding  the  method 
of  color  photography  used  here. 

A.    THE  PHOTOGRAPHIC  PROCESS* 

A  conventional  photographic  film  comprises  a  structural 
substrate  on  which  is  placed  a  colloidal  suspension  of 
crystals  of  one  or  more  of  the  silver  halides:  silver 
bromide,  chloride  or  iodide.  Each  crystal  comprises  silver 
and  halide  ions  in  a  predominantly  regular  crystalline 
structure.  Within  the  structure,  however,  there  are  a  few 
positive  (Ag*)  silver  ions  which  occupy  interstitial 
positions.  When  a  a  packet  of  electromagnetic  radiation  (a 
photon)  of  sufficient^  energy  impinges  on  a  crystal,  one  or 
more  electrons  are  broken  loose  from  the  halide  ions.  This 
results  in  a  modification  of  the  state  of  the  crystal.  The 
few  free  electrons  diffuse  through  the  crystal  until  they 
come  in  contact  with  some  of  the  interstitial  silver  ions. 
The  positive  silver  ions  capture  the  electrons  and  become 
metallic  silver  atoms.  A  crystal  with  one  or  more  metallic 
silver  atoms  is  termed  a  "sensitized"  crystal.  On  a 
macroscopic  scale,  the  distribution  of  these  sensitized 
crystals  provides  an  image  of  the  incoming  light  pattern. 
Since  the  crystals  have  not  changed  entirely  to  metallic 
silver  the  image  is  not  normally  visible,  and  is  termed  a 
"latent"  image.    When  exposed   to   a   developer,   however, 


♦Engel,   Charles   E. ,   Photography   for   the  Scientistx 
Academic  Press,  1968  ,  p.2.ffT 
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the  atoms  cf  silver  catalyze  the  reaction  of  the  developer 
with  the  remainder  of  the  crystal,  causing  it  to  convert 
entirely  to  metallic  silver,  and  thus  forming  what  is 
commonly  known  as  the  negative.  Subsequently,  the  negative 
is  exposed  first  to  chemicals  which  neutralize  the  developer 
(a  "stop"  hath)  and  then  to  those  which  remove  the 
unsensitized  crystals  from  the  negative  (a  "fixer"  bath) . 

The  chemical  reactions  described  all  reguire  energy.  In 
particular  the  crucial  event,  that  of  the  halide  ion  giving 
up  an  electron,  as  in  the  bromide  ion, 

Br j — *Br  +  e- 

energy 

only  occurs  when  the  applied  energy  is  sufficient  to  break 
the  molecular  bond.  The  energy  required  is  proportional  to 
the  oxidation  potential  of  the  element  and  is  found  to  be 
approximately  10-19  joule.5 

B.    FREQUENCY  AND  ENERGY  IN  PHOTOGRAPHY 

The  energy  required  to  break  electrons  free  from  a 
bromide  ion  must  be  provided  by  the  photon  striking  the 
crystal.  The  energy  of  a  photon  is  known  from  quantum 
mechanics  to  have  energy  proportional  to  its  frequency 

E  =  hf 

where  h  is  Planck's  constant,  6.626  x  10-34  joule-sec.  and 
f  is  the  frequency  of  the  radiation.6  The  lowest  frequency 
radiation  which  is  theoretically  capable  of  triggering  the 
process  is  then  found  to  be  approximately 

E  =  f  =  10-19  =  1.60  x  10i*  Hz 


SMasterton,   W.   L.   and   Slowinski,   E.   J.,   Chemical 
Principles,  W.  B.  Saunders,  1973,  p. 255. 

6Beiser,    Arthur,    Perspectives   of   Modern   Physics,. 
McGraw-Hill,  1969,  p.  56. 
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or  160  thousand  gigahertz,  which  corresponds  to  a  wavelength 
of  2  microns  (2  x  10~6  meters) .  (The  maximum  wavelength 
that  can  be  imaged  by  available  Kodak  films  is  reported  to 
be  1.15  microns.7)  By  contrast,  the  radiation  associated 
with  the  radio  spectrum  is  of  the  order  of  ten  gigahertz  or 
less,  four  orders  of  magnitude  lower  than  the  necessary 
frequency  level.  Thus,  it  is  not  likely  that  direct 
irradiation  by  RF  energy  shall  trigger  the  photographic 
process. 

B.    THE  COLOR  PROCESS9 

Imaging  fields  is  possible  however,  by  controlling  the 
events  which  occur  subseguent  to  the  sensitization  of  the 
film.  In  particular,  the  Polaroid  Corporation  color  process 
("Polacolor")  uses  a  diffusion  technique  which  is  sensitive 
to  temperature.  If  the  RF  energy  field  can  be  transduced 
into  a  temperature  field,  modulation  of  the  diffusion 
process  is  possible,  resulting  in  an  image  of  the  field 
being  formed. 

The  principle  of  crystal  sensitization  in  color  film  is 
the  same  as  described  above  for  black  and  white  film.  The 
difference  between  color  and  black  and  white  film  is  that 
color  film  actually  has  three  emulsions,  while  black  and 
white  film  has  only  one.  Each  of  the  three  layers  is 
sensitized  by  a  different  primary  color,  the  undssired 
colors  having  been  blocked  by  filters  built  into  the 
emulsion.  In  addition  to  the  silver  halide  crystals  and  the 
color  filter,  each  layer  has  a  chemical  (a  "dye-coupler") 
which  reacts  with  the  sensitized  crystals  and  the  developer 
to  create  a  dye  when  the  film  is  developed.  It  is  at  this 
point  where  the  photographic  process  can  be  used  to  provide 
imaging  of  electromagnetic  energy.  Unlike  other  processes, 
the   Polacolor   method  causes  dyes  to  be  formed  in  all  three 


'Eastman  Kodak  Company,   Kodak   Publication   No.   N-17, 
Kodak  Infrared  Films,  1972,  p.  2. 


Engel,  p.  45 ,  ff. 
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layers  of  the  film  independently  of  the  amount  of 
sensitization  of  the  emulsion.  Subsequently,  sensitized 
crystals  in  each  of  the  three  emulsions  cause  some  of  their 
respective  dyes  to  become  anchored  in  the  emulsion.9  This 
anchoring,  or  leaving  behind  of  selected  dyes  is  termed 
mordanting.  How  much  of  each  dye  is  mordanted  (rather  than 
the  amount  created)  is  therefore  proportional  to  the  amount 
of  sensitization  to  which  the  film  has  been  exposed.  The 
unanchored  dye  molecules  then  diffuse  upward  through  the 
multilayered  emulsion  until  reaching  the  substrate  for  the 
print,  where  they  become  fixed.  It  is  this  diffusion 
process  which  is  used  in  dye-transfer  modulation  to  obtain 
an  image. 


«Thirtle,  John  R.,  "Chemistry  of  Color  Process,"  SPSE 
Handbook  of  Photographic  Science  and  Engineering,  Wiley, 
TO73-p7567-56  87 
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III.   THE  PROCESS  OF  DYE  TRANSFER  MODULATION  (DTM) 

The  Dye  Transfer  Modulation  (DTM)  Process  is  primarily 
dependent  on  the  Polacolor  diffusion  process.  The  technique 
comprises  four  sub-r  processes. 

A.    BASIC  PRINCIPLE  OF  DTM 

It  can  be  shown10  that  the  rate  of  a  chemical   reaction 

is  temperature  dependent  by  the  expression 

R        T   -  T 
log   — 2.  =  k  — 2 i- 


10  r      ITT 

1  12 

where   K   is  a   constant,   T    is  the   temperature   of   the 

1  (i) 

reaction   and   R  is   the   rate   of    reaction   at   the   correspond- 

(i) 

ing  temperature.  It  can  further  be  shown11  that  the  rate  of 

diffusion   of   a   substance   through  a  medium  is  temperature 
dependent  by  the  expression 

T   -  T 


j  =  K 


2     T  T 
1  2 


where   j   is   the  change  in  diffusion  flow  for  the  indicated 
change  in  temperature  T  and  K   is  another  constant.   It  can 

2 

be  seen  that  for  both  rates  the  effect  of  a  given  change   in 

temperature   at   a   low  temperature  is  significantly  greater 

than   at   a   higher   temperature.   For   example,    at    room 

temperature   a   five   degree   change  in  reaction  temperature 

gives  the  ratio 

T   -  T 
2     i    =  5       5.46  x  10-s 

T  T      300*305 

1   2 

while   at   the   temperature  of  dry  ice,  approximately  195°X, 

that  same  five  degree  change  in  temperature  yields 

T   -  T 

__2 I  5     =  1.28  x  10-* 

T  T      195*200 

1   2 


iOMasterton   and   Slowinski,    p. 380-382. 
^Encyclopaedia    Britannica,     1968.    v. 7,    p. 424 
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2.35  larger  than  the  effect  at  room  temperature.  It  is 
clear  that  temperature  variations  can  be  used  to  modulate 
the  reaction  rate  in  the  formation  of  dyes  and,  further,  to 
modulate  the  rate  at  which  the  dyes  in  the  unsensitized 
layers  diffuse  to  the  print  substrate.  It  is  also  clear 
that  a  given  change  in  temperature  shall  have  a  greater 
effect  on  these  rates  if  the  ambient  temperature  is  low. 

The  basic  principle  of  the  Dye  Transfer  Modulation 
Process  is  to  transduce  an  energy  field,  particularly  an 
electro-magnetic  field,  into  a  thermal  field,  then  use  the 
thermal  field  to  modulate  the  rates  of  reaction  and 
diffusion  in  the  Polacolor  process.  The  term  "Dye  Transfer 
Modulation"  reflects  the  fact  that  the  end  result  is  to  vary 
the  amount  of  dye  which  transfers  from  the  film  negative  to 
the  print. 

B.    FILM  SENSITIZATION 

The  schematic  drawing  of  Figure  1  shows  the  cross'- 
section  of  a  Polaroid  Film  Pack.  Figures  2  and  3  show  the 
actual  film  packet  and  its  components.  Three  emulsion  layers 
are  placed  on  the  negative  substrate,  sensitive  to  red, 
green,  and  blue  light,  in  order  upward  from  the  substrate. 
Their  developers,  when  activated,  produce  cyan,  magenta,  and 
yellow  dyes. 12  Thus  the  layer  farthest  from  the  surface  of 
the  negative  is  the  layer  sensitive  to  red  light,  which 
produces  cyan  dye.  Controllability  of  the  diffusion  process 
is  maximized  by  allowing  only  this  deepest  dye  to  escape  the 
negative.  This  is  accomplished,  in  principle,  by  exposing 
the  film  thoroughly  to  blue  and  green  light.  Such  exposure 
will  cause  the  yellow  and  magenta  layers  to  be  mordanted  in 
the  negative,  leaving  only  the  cyan,  red-sensitive  layer 
free  to  diffuse  to  the  surface  of  the  negative  and  thence 
onto  the  print  surface. 

This  selective  pre-sensitization  of  the  film  is 
effected   by   exposing   the   film   in  a  box  camera  and  using 


X2SPS2  Handbook  .. . ,  p.  567 
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Figure  2 
Polaroid  Film  Pack 


Figure  3 
Exploded  View 
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filters  to  allow  only  the  desired  wavelength  ranges  of  blue 
and  green  light  to  reach  the  film.  The  sensitization  setup 
is  shown  in  Figure  4.  Practical  aspects  and  parameters  of 
this  step  in  the  process  are  discussed  below. 

C.  TEMPERATURE  BIASING 

The  film  is  then  placed  in  a  closed  styrofoam  container 
(shown  in  Figure  5)  at  a  position  approximately  one 
centimeter  away  from  a  layer  of  dry  ice  and  allowed  to  cool 
to  a  selected  temperature. 

D.  DEVELOPER  ACTIVATION 

When  properly  cooled  the  developers  in  the  film  are 
activated  by  passing  the  film  pack  through  a  Polaroid  4x5 
Land  film  holder  model  500  (figures  6  and  7) .  The  holder  has 
rollers  which  break  a  pod  containing  an  alkaline  activation 
chemical  (shown  in  figure  3)  and  spread  the  activator 
between  print  and  emulsion.  The  activator  diffuses  into  the 
emulsion  and  activates  the  developer. 

E.  TRANSFER  MODULATION 

The  film  pack  is  then  placed  in  an  electromagnetic 
field.  Here,  current  is  induced  in  the  emulsion  causing 
localized  heating  proportional  to  the  intensity  of  the 
electric  field.13  This  heating  then  speeds  the  diffusion  of 
the  cyan  dye  molecules  to  the  print  substrate  producing  an 
image  of  cyan  dye  in  the  areas  of  greater  field  intensity. 


13Keigo  Iizuka,  p.  3. 
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Figure  4 
Film  Sensitization  Setup 


Figure  5 
Film  Cooling 
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Figure  6 
Polaroid  Film  Holder 


Figure  7 
Film  Activation 
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IV.  PARAMETERS  OF  THE  PROCESS 
The  thesis  objective  was  to  optimize  the  dynamic  range 
of  the  process.  The  dynamic  range  is  primarily  determined 
by  the  film  characteristics  and  by  the  various  process 
parameters  which  affect  those  characteristics.  This  chapter 
identifies  and  qualitatively  describes  some  of  these 
parameters. 

A.  FILM  DENSITY 

A  film  characteristic  of  fundamental  value  here  is  the 
concept  of  photographic  density  (D) .  Density  is  a  measure 
of  the  darkness  of  a  film  image.  It  is  quantitatively 
defined  in  sensitometric  terms.1*  In  the  Polacolor  process, 
how  dark  the  image  gets  depends  on  the  amount  of  coloring 
matter  transferred  to  the  print.  As  described  above,  the 
film  is  subjected  to  a  radiation  field  after  the  dye 
transfer  process  has  been  activated.  Initially,  the  print 
surface  is  devoid  of  any  color.  As  the  film  packet  is  warmed 
by  the  field  and  by  the  ambient  temperature  the  dyes  in  the 
cyan  layer  (which  has  not  been  exposed  to  light)  begin  to 
transfer  to  the  print  by  diffusion  through  the  activator 
chemical  which  is  between  them.  Obviously,  if  the  transfer 
process  is  allowed  to  continue  for  a  longer  time  more  dye 
molecules  can  reach  the  print. 

B.  CONTRAST 

When  measured  experimentally,  density  becomes  a 
space-averaged  value  and  does  not  provide  any  information  on 
the  relative  distribution  of  light  and  dark  material  or  of 
color  in  the  print.  Yet  the  resolution  and  clarity  of  the 
image  are  apparently  highly  dependent  on  both  these  factors. 
The  quantity  contrast  provides  a  measure  of  the  perceived 
distribution  of  light  and  dark  regions  in  a  print.  It  is 
essentially  a   subjectively   determined   quantity,   although 

l*Appendix  A  p.  87. 
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methods  exist  to  represent  the  contrast  characteristic.  The 
industry  uses  gamma  (  "f  )  where  T  is  the  slope  of  the 
straight  line  portion  of  the  characteristic 
density.  yergus_exposure  ("D  -  Log  H")  curve.15  Gamma 
represents  the  maximum  amount  of  contrast  for  a  given  curve. 
The  industry  also  uses  Contrast  Index  (CI)  as  a  measure  of 
contrast. 16 

C.    APPLIED  ENERGY 

An  input  parameter  of  prime  significance  is  that  of 
applied  energy.  As  seen  from  photographic  science,  density 
is  highly  dependent  on  the  applied  energy.17  Since  energy 
is  simply  power  times  time  in  constant  power  situations, 

E  =  P«T 

it  is  concluded  that  one  may  relate  density  to  power  being 
expended.  The  quantity  generally  used  is  the  rate  of 
energy  flow  per  unit  area,  or  p_ower  density,  measured  in 
watts  per  square  centimeter.  This  quantity  was  found  to  be 
significantly  dependent  on  the  interaction  of  several 
contributory  factors: 

1 •   Power^Source 

The  power  applied  to  the  DTM  process  depends  most 
importantly  on  the  sources  available.  Selection  must  be 
made  among  these  on  the  basis  of  still  more  factors.  In 
general  the  sources  available  were  limited  to  approximately 
300  mw/cm2.  This  limitation  is  due  primarily  to  the  factor 
discussed  next. 

2-   Source  Frequency 

The  operating  frequency  is  bounded  below  at 
approximately  1.5  gigahertz.  This  bound  is  set  by  the 
dimensions   of   the   film  and  of  the  target  structures  whose 

15  See  Appendix  A,  p.  87. 

i&Kodaic  publication  #P-315  "KODAK  Plates  and  Films  for 
Scientific  Photography",  Eastman  Kodak  Company,  1973, p.  3. 

17See  Appendix  A,  p. 87. 
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field  characteristics  are  to  be  investigated.  An  upper 
limit  is  set  by  the  ability  to  visually  differentiate 
features  in  a  pattern.  A  selection  was  possible  among 
sources  which  have  output  freguencies  between  .2  -  10 
gigahertz.  Lower  freguency  devices  can  deliver  larger  power 
levels,  but  their  radiating  devices  and  target  elements  are 
generally  larger  than  those  of  the  higher  freguency 
eguipraent,  so  that,  with  the  available  sources,  the  area  to 
be  irradiated  varied  directly  with  the  power  available, 
resulting  in  an  essentially  constant  availability  of  power 
density.  Conseguently  freguency  selection  was  made  only  on 
the  basis  of  the  operating  freguency  bounds.  The  freguencies 
ultimately  selected  were  approximately  1.5  and  8   gigahertz. 

3«   Sou r ce_and_ Tar get _Re so nances 

A  second  technologically  related  factor  is  that  of 
resonance.  The  limitations  on  available  power  density 
reguired  that  experiments  be  conducted  in  the  near  fields  of 
the  radiating  devices.  This  situation  results  in  strong 
interaction  of  mutual  impedances,  severely  complicating  the 
task  of  matching  the  source  impedance  to  its  load.  The  time 
response  of  this  near-field  system  has  time  constants  as 
shown  by  equation  V.  C.  23  on  page  60.  The  time  constants  in 
these  experiments  were  sufficiently  high  in  some  cases  to 
cause  the  impedances  involved  apparently  to  take  on  a 
time-varying  aspect.  That  further  complicated  the  task  of 
maximizing  the  power  delivered  to  the  film. 

(| .   Source  Radiating_Devices 

The  different  radiating  devices  reguired  for  the 
various  sources  result  in  different  field  distributions 
(and,  conseguently,  differing  power  densities)  being  applied 
to  the  target  elements.  Ultimately  two  devices  were 
selected  which  maximized  power  density  available  at  the 
selected  freguencies.  A  guarter- wavelength  monopole  over  a 
ground  plane  with  a  corner  reflector  was  used  for  the  lower 
freguency  work,  while  an  open  ended  waveguide  was  used  at 
the  upper  freguency. 
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5-   Target  Element  Placement 

A  final  power-related  aspect  is  that  of  target 
element  placement,  which  affects  interaction  among  the 
various  elements  of  the  near  field  scenario  and  the 
consequent  resultant  load  impedance  presented  to  the  source. 

D.    PHOTOGRAPHIC  PARAMETERS 

The  photographic  parameters  are  those  which  affect  the 
print  through  mechanisms  primarily  related  to  photographic 
science. 

1 .   Color Temperature   of Liqh  t_ Source and Zii£ 

Sensitivity 

Appendix  A  provides  an  overview  of  quantitative 
photographic  terms  and  concepts  which  shall  be  used  in  the 
discussion  which  is  to  follow.  Figure  8  provides  a 
representative18  sensitivity  plot  for  Polaroid  type  58  film 
and  indicates  peak  sensitivities  to  blue,  green,  and  red 
light  (of  the  yellow,  magenta,  and  cyan  layers  respectively) 

of 

Blue  350nm. 

Green  350nra.,  572nm. 

Red  645nm. 

It  is  particularly  interesting  to  note  that  the  green- 
sensitive  layer  has  two  sensitivity  spectra,  one  which  is  a 
proportional  mapping  of  the  same  spectrum  to  which  the  blue 
layer  is  reported  to  be  sensitive.  It  is  stressed  that  the 
dyes  do  not  function  completely  independently: 

Although   each  dye  absorbs  in  one  principal  spectral 
reqion...,   it  also  has  unwanted  density  in  the 
other  spectral  regions.   These  unwanted  densities 
desaturate  and  darken  colors  in  a  photograph,  thus 
distorting  color  reproduction.19 

This   duality   in   spectral   sensitivity,    and   the 

relatively   large  separation  between  the  blue-green  spectrum 


18Polaroid  Corporation  letter  received  March  1974  and 
appended  as  Appendix  B.  Note  the  caveat  that  these  curves 
are  not  to  be  considered  as  average  or  a  product  specifica- 
tion. Data  confirmed  as  being  equally  applicable  to  type  58 
as  to  type  108  by  telephone  conversation  21  June  1974  with 
Polaroid  Corporation  Technical  Assistance  Division 
representative  Mr.  Brooks  Corl,who  reported  that  the  type  58 
ana  108  films  use  the  same  identical  emulsion. 

19SPSE  Handbook. . . ,  p.  452-453. 
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Figure    8 
Sensitivity    of    Polaroid   Type    Film 
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Figure  9 
Passbands  of  Selected  Wratten  Filters 
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and  the  red-sensitive  region  tend  to  indicate  use  of  the 
red-sensitive  layer  for  the  dye  image,  since  the  blue  and 
green-sensitive  layers  should  be  able  to  be  mordanted  in  the 
negative  without  excessively  activating  the  red-sensitive 
layer. 

The  color  temperature  corresponding  to  the  color  to  be 
used  for  sensitizing  the  film  are20 

Color     Zavelengthlaajl     Color  Temperature  (°K) 

Blue  350  8279 

Green  350  8279 

572  5066 

These  colors  are  provided  approximately  by   Wratten   filters 

18A   (at   360   nm)  and  99  (at  550  nm) .    See  Figure  9.    The 

use  of  filters  such  as  these,   whose   pass-band   half-widths 

are  45  and  35  nm.,   respectively,21  minimizes  the  problem  of 

trying  to  match  color  balance  of   the   film   as   long   as  a 

constant  intensity  of  sensitizing  light  is  used.   (Increased 

exposure  at  any  particular  wavelength   gives   less   negative 

density.22     Thus,   sufficient   increase   in   exposure   may 

override  the  attenuation  characteristic  of  the  filter.) 

If   large   variations   in  exposure  are  to  be  used,  then 

the  light  source  must  be  matched  to   the   color   temperature 

balance   of  the  film  using  a  color  conversion  filter  with  an 

appropriate  Mired-Shift-Value   (MSV).23    Polaroid   film   is 

color   balanced   for   5500°K.2*    For   use   with  4800°K  blue 

photoflood  lamps  a  color  conversion  filter  with   an   MSV   of 

-26.5   is   indicated,   which   corresponds  approximately  to  a 

Wratten  78C  (whose  MSV  is  -24). 25 


20Appendix  A,  p.  88. 

2lSPSE  Handbook ,  p.  308-309. 

22c.f.  Appendix  A,  p.  87. 

23Ioidi,  p.  9  2. 

2*Polaroid   Corporation   publication   PX   850,  December 
1972,  p.  1. 

2  5SPSE  Handbook ,  p. 315. 


28 


In  summary,  the  cyan  dye  layer  is  indicated  to  be  the 
most  useful  layer  for  the  modulation  process.  The  yellow 
and  magenta  layers  should  be  mordanted  in  the  negative  by 
exposure  cf  the  film  to  blue  light  of  350nra.  and  green  light 
of  572nm.  approximately  derived  by  exposing  the  film  to  a 
U800°K  light  source  first  through  a  Wratten  18A  blue  filter, 
then  a  Wratten  99  green  filter.  If  various  exposure 
settings  are  anticipated  the  light  source  should  be 
compensated  to  5500°K  by  concurrent  additional  use  of  a 
Wratten  78C  color  conversion  filter. 

2 .   Choice r of ^Co lor _Exgosure_Ti mes_f or ^Sensitization 

The  film  negative  has  been  described  as  comprising 
three  basic  dye  layers.  The  pre-exposure  of  the  film  is 
seen  to  be  a  method  of  selectively  mordanting  dye  layers  in 
the  negative  to  allow  only  the  deepest  (cyan-colored)  layer 
to  diffuse  up  to  the  print.  It  is  not  considered  likely 
that  in  practice  the  two  higher  layers  are  completely 
mordanted.  Rather,  it  is  apparent  that  to  some  extent  the 
yellow  mixes  with  cyan  to  provide  a  green  tone,  while 
magenta  mixes  to  provide  blue.  The  relative  proportions  of 
these  can  be  controlled  by  varying  the  amount  of  light  each 
receives  during  pre-exposure.  Increasing  blue  light  results 
in  less  yellow  being  transferred.  Increasing  green  light 
results  in  less  magenta  being  transferred. 

As  a  consequence,  rather  than  a  single  D-log  H 
curve,  there  are  two  density  curves  in  this  process;  density 
versus  blue  exposure  for  selected  values  of  green  exposure, 
and  density  versus  green  exposure  for  selected  blue 
exposure.  It  is  possible  to  select  values  of  blue  and 
green  pre-exposure  which  set  an  operating  point  on  the  three 
dimensional  contour  which  provides  the  largest  value  of 
<ja  m  ma . 

These  curves  are  found  also  to  vary  with  respect  to 
the  initial  temperature  selected,  as  discussed  below. 
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3.   Developer  Activation 

The  film  packet  must  be  pulled  through  the  film 
holder  at  a  rapid  and  uniform  rate.  Any  hesitation  while 
pulling  the  packet  through  the  rollers  results  in  uneven 
distribution  of  the  dye-developer  activator  with  a  resultant 
ridging  in  the  final  image. 

E.    THERMAL  PARAMETERS 

1 .  Temperature  Biasing 

The  film  packet  must  be  sufficiently  cooled  to 
ensure  that  the  only  modulation  of  diffusion  and  development 
which  occurs  is  due  to  the  localized  heating  by  the  incident 
energy  field.  If  chilled  excessively  reactions  cannot 
occur.  Further  cooling  of  the  reagent  pod  will  cause  it  to 
solidify  and  thus  prevent  activation  of  the  developers. 
Conversely,  if  the  packet  is  not  sufficiently  cooled  the 
development  and  dye  transfer  processes  will  proceed  at  a 
rate  too  high  to  allow  significant  modulation  by  the 
incident  field.  The  effect  of  these  varying  factors  is  seen 
primarily  as  a  change  in  print  density.  It  can  be 
understood  that  if  the  film  is  cooled  for  different  amounts 
of  time  or  with  different  materials,  its  temperature,  when 
initially  inserted  in  the  radiation  field,  shall  be 
different.  For  a  lower  initial  temperature  a  film  should  be 
expected  to  require  more  heating  (longer  radiation  time  or 
higher  incident  energy  density)  to  acquire  a  particular 
density  reading.  It  may  be  then  concluded  that  the  print 
density  is  dependent  on  initial  temperature  as  well  as  on 
applied  energy. 

2.  Cooling  .Procedure _ Para meters 

Several  aspects  of  the  cooling  technique  were 
investigated,  primarily  as  a  means  of  calibrating  the 
temperature  of  the  film. 

a.   Choice  of  Coolant 

Two  methods  were  used  to  cool  the  film.   Liquid 
nitrogen   was   used   for   very   low   temperature   work    (at 
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temperatures   below   -30°C)  ,   and   dry   ice  was   used   for 
temperatures  between  -30°  and  20°  C. 

b.  Cooling  Device  Temperature  Variations 

These  were  found  to  be  essentially   exponential 
responses  in  time,  and  highly  predictable. 

c.  Activator  Chemical  Freezing  Point 

This    value    was   investigated   in   order   to 
determine  a  lower  bound  on  the  temperature  for   pre-chilling 
the  film. 
F.    DYNAMIC  RANGE 

For  a  given  setting  of  other  process  parameters,  the 
amount  of  density  (D)  appears  to  be  proportional  to  the 
amount  of  energy  (S)  present  in  the  fields  being  measured. 
Therefore  one  measure  of  dynamic  range  might  be  the  range  of 
energies  over  which  noticeable  changes  in  density  occur. 
Referring  to  the  D  -  log  H  curve26  one  can  see  that  dynamic 
range  can  be  measured  by  the  range  of  energies  over  which 
there  is  a  slope  to  the  curve. 
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Typical  D  -  Log  H  Curve 


26Appendix  A  p.  87. 
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Alternatively,  it  can  be  seen  that  the  dynamic  range  is 
also  affected  by  the  total  amount  of  change  in  density.  A 
large  range  of  energies  is  of  little  use  if  there  is  only  a 
slight  change  in  the  density  over  that  range.  Therefore  it 
is  concluded  that  the  range  of  values  over  which  the  film 
density  varies  for  different  energy  settings  is  also  an 
important  measure  of  the  dynamic  range  of  the  process. 

A  performance  measure  of  the  dynamic  range  of  the 
process,  which  takes  into  consideration  ranges  of  density 
and  energy,  is  a  quantity  defined  as  sigma  (c): 


o-*   |j*  M>SE  «  AD  AE 


It  is  considered  desirable  to  find  the  settings  of  other 
process  parameters  which  result  in  an  optimization  of  sigma 
over  its  range  of  values. 
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v-    EXPERIMENTAL  METHODS  AND  MODELS 

Several  characteristics  and  parameters  affecting 
optimization  of  the  process  have  been  identified.  The 
quantification  of  these  factors  was  achieved  through  various 
experiments.  The  following  is  a  description  of  those 
experiments  and  of  the  models  used  to  develop  them. 

A.    DETERMINATION  OF  FILM  DENSITY 

1 •   Film, Density^Standard 

Film  density  was  deemed  a  fundamental 
characteristic.  An  ability  to  measure  it  was  considered 
essential  for  other  experiments.  Facilities  were  not 
available  for  direct  density  measurements  in  small  areas  of 
the  film.  Therefore  a  set  of  prints  was  developed  as  a 
comparison  standard  for  selected  combinations  of  blue  and 
green  sensitization.  Other  photographs  could  then  be 
compared  visually  with  the  standards  in  order  to  determine 
density  values  in  different  regions. 

The  standards  were  developed  as  follows: 

a.  A  particular  ratio  of  blue  and  green 
sensitizing  exposure  was  selected;  e.g.  2:1. 

b.  Successive  exposures  were  made  of  a  uniform 
surface  (white  cardboard)  with  a  constant  light  source 
(4800°K  at  225  foot-candles.)  Each  exposure  was  made  with 
the  selected  ratio  but  with  increasing  fractions  of  seconds 
of  exposure  for  each  color:  e.g.:  1sec.  blue    1sec.    green; 

s  10 

2sec.  blue  lsec.    green;  1sec.  blue   lsec.  green;  etc. 

5  5  2 

c.  Each  exposure  was  developed  at  a  uniform 
temperature  (23°  C)  for  60  seconds. 

d.  A  method  was  devised  to  provide  the  opacity  (0) 
reading: 

Incident  Intensity  (I  ) 
0  = i 

Reflected  Intensity  (I    ) 

ref 
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from  which  the  density  (D)  was  computed 

D  =  log   (0) 

10 

The  method  used   fluorescent   light   on   a   white   cardboard 
surface    to   reflect   uniform   light   on   the   print.   The 

light  (I  )  incident  on  the  print  from  the  cardboard  and   the 

i 
light  (I    )  reflected  from  the  print  were  measured   with   a 

ref 
light   meter27   in   foot-candles.     Figure   11   shows   the 

measurement   method.    Densities   (D)  were  obtained  for  each 

exposure   for  three   different  values  of  I   (38,  50,  and   60 

i 
foot-candles) .   These  values  were  then  averaged   to   give   a 

value  of  average  density  (D) . 

e.   A    different   ratio   of   blue   and   green   was 

selected  and  the  process  repeated. 


Figure  11 
Density  Measurement  Method 

2 .   Variation of Density with Preexfio sur e Light 

Intensity  ..and  Color_Balance 

The  combined  effect  of  intensity  and  color  balance 
on  film  density  was  obtained  by  plotting  the  various  values 
of  Density  obtained  above  in  a  three  dimensional  display  of 
D  versus  Log  H  for  blue  and  green  preexposure. 

27General  Electric  Model  8DW58YI 
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B.    DETERMINATION  OF  DYNAMIC  RANGE 

Another  fundamental  characteristic  of  the  process  to  be 
measured  was  the  Dynamic  Rar^e,  the  breadth  of  which  is  to 
be  characterized  by  the  quantity  cr  (siqma) : 

where  M)  is  the  range  of  print  densities  resulting  from  a 
range  of  applied  energy  AE.  The  value  of  this  performance 
measure  was  to  be  found  for  various  ratios  of  blue  and  green 
light  sensitization,  for  different  levels  of  total  applied 
energy  and  for  different  settings  of  other  process 
parameters.   These  quantities  were  measured  as  follows: 

1.   A   particular   ratio  of  blue  to  green  sensitization 

exposure  was  selected. 

2  The  film  »as  pre-chilled  to  a  specifiad 
temperature,  activated  and  placed  at  the  open  end  of  a 
waveguide.  The  power  supply  was  turned  on,  irradiating  tne 
film  for  a  specified  interval  of   time   during   which   power 

levels  were  recorded. 

3  At  the  end  of  the  radiation  period  the  power 
supply  was  turned  off  and  the  print  was  separated  from  the 
negative  to  halt  the  transfer  process. 

<,  If  an  image  was  visible  in  the  region  of  the  wave- 
guide aperture  then  the  power  was  reduced  using  a 
calibrated  attenuator.  The  process  was  repeated  until  there 
was  a  just  barely  noticeable  change  in  density.  St  that 
point  the  highest  energy  level  Emax  of  the  energy  range  was 
considered  to  have  been  reached. 

5  The  width  of  the  image  was  measured  between  points 
at  which  an  image  was  gust  barely  visible.  This  width 
could  be  converted  to  a  power  difference  with  respect  to  the 
maximum  power.-  At  that  point  it  was  considered  that  the 
lowest   energy   level  Emin  had  been  reached  for  the  selected 

28  (see  paragraph  V  C  below). 
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ratio.     Therefore   a   new   ratio   was   selected   and  the 
experiment  repeated. 

6.  When  all  desired  ratios  had  been  tested  the  recorded 
power  levels  were  adjusted  for  factors  such  as  reflections29 
to  find  the  actual  power  delivered  to  the  film.  The 
quantity  AE  was  then  found  to  be  the  difference  between  the 
adjusted  highest  power  level  Emax1  and  the  adjusted  lowest 
power  level  Emin': 

&E  =  Emax'  -  Emin* 

7.  The  prints  at  Emax'  and  Emin1  were  measured  for 
density  by  comparison  with  the  density  standards.  The 
difference  in  density  at  the  two  respective  energy  levels 
was  taken  to  be  aD . 

C.    DETERMINATION  OF  APPLIED  ENERGY 

Since  the  object  of  this  work  was  to  image  applied 
fields,  it  was  necessary  to  establish  a  means  of  measuring 
the  power  applied  to  the  film  by  the  selected  sources. 
There  are  two  aspects  to  this  problem.  One  is  measuring  the 
power  applied  to  the  film  placed  at  the  end  of  the  waveguide 
for  dynamic  range  measurements.  The  other  is  measuring  the 
power  applied  to  the  film  placed  in  the  field  of  an  array  of 
dipoles.  In  either  case  energy  E  is  taken  to  be  the  product 
of  time  averaged  power  (P)  and  the  radiation  period  (T) . 

E  =  P«T 

In  both  cases  there  were  system  time  constants  which 
resulted  in  a  slow  change  in  power  delivered.  These 
constants  are  considered  to  be  due  primarily  to  impedance 
characteristics  of  the  film30  since  pains  were  taken  to 
maintain  the  power  supply  frequency  and  power  output 
constant.  Because  of  these  time  constants  the  steady  state 
power  levels  were  taken  as  the  values  to  be   recorded.    The 

3° (see  paragraph  V.  C.  3  below) . 
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following    derivations   therefore   assume   steady   state 
conditions. 

The    schematic  of  the  power  measurement  setup  is  shown 
in  Figure  12. 

POWER 
REFLECTED    FORWARD 


POWER 
METER 

& 
MOUNT 


POWER 

METER 

S 

MOUNT 

V 


calibrated  attenuators 


rDirectional\ 
i  Coupler 


POWER 
SUPPLY 


Figure  12 
Power  Measurement  Setup 

In  both  cases  power  was   measured   through   directional 

couplers,  thermistor  mounts  and  power  meters.   The  equipment 

used  were: 

Hewlett  Packard: 

Model  U30C  Power  Meter 

with 

Model  U77B  Thermistor  Mount 

766D  Directional  coupler 
Assorted  calibrated  attenuators 

1  ■   PQ,1:e£!°iJ13.£i9Q.  Q£ _Po wer_ Deli vered_bY_ Waveguide 
a.   Field  Equations 

For  a  rectangular  waveguide  (Figure  13)  driven 
by  a  sinusoidal  source  of  such  a  frequency  (w)  as  to  induce 
TE    mode   operation   above   cutoff,   the   solution   to  the 

10 

waveguide   equations  gives  the  following  expressions  for  the 
X#  Y,  and  Z  components  of  the  electric  (E)  and  magnetic   (H) 
fields: 3i 


31Jordan,  E.C.    and  Balraain,  K.G.,  Electromagnetic  Waves 
and  Radiating*  Systems,  2nd  Ed,  Prentice-HaTT7~T96H,  pT-?^. 


37 


iwt 
^v    -    r>  Hx   =  rtrHo     sin  irx   e  x 

Ex    "    u  a"K~*  a"* 

jotA  _* 

Ev  =  -j  ttfiTtH*    sin  -rrx   e        y  Hy   -    0 

5^        a  „     .  j-t . 

EZ    =    0  Hz    =    He    cos  TLx    e  z 

a 

where    for   TE        mode    above   cutoff 
1  o 

Kc2    =   (in2    =    T2    +  <o2yU£ 

a   =  the  X  dimension  of  the  waveguide  :-j 
IT  =  an  arbitrary  complex  constant  pjej 

al  and  £  are  respectively  the  permeability  and  permittivity 
constants,  and  2,  y,  z  are  unit  vectors  in  rectangular 
coordinates.  When  expressed  in  phasor  form  these  simplify  to 

fx  =  0  H x  =  j^J.  sin(-|x)x 

Ey  =  -jo^ej  sin(Tmy    Hy  =  0  (V-  C.  1) 

fz    =    0  Hz  =  lHol  cos(|xjz 

b.   Time  Averaged  Power  Density 

The   time   averaged   power   density  is  known  to 

be32  a  vector  guantity 

p  =  i  Re  {Ex  H*}  (v-  c-  2) 

2 

When  the  cross  product  E  x  H*  is  taken,  the  real  component 
is  simply  E^hSz  so  that  time  averaged  power  flow  (P(x))  in 
watts  per  sguare  meter  is 

P(X)    =   i(EyHx)    Z 

=  ^afj^li  sin2(|x)z     (V.  C.  3) 

Figure  13  shows  this  distribution  as  seen  across  the  face  of 
the  waveguide.  It  shows  a  flow  of  power  down  the  waveguide 
in  the  z  direction,  uniformly  distributed  in  the  y-direction 
but   tapering  off   sinusoidally   as   it  moves  away  from  the 

3  2ibicU,  p.  171. 
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center  in  the  x-direction.  Figure  14  is  a 
photograph  of  the  waveguide  cross  sec  .ion  showing  that 
distribution  in  fact. 


n&w«c 


*=a  "X 


Figure  13  Waveguide  Dimensions  and  Power  Distribution 


Figure  14 
Photograph  of  Power  Distribution  across  Face  of  Waveguide 

Clearly,  there  is  in  fact  a  power  distribution  across  the 
waveguide.  The  problem  is  how  to  relate  those  levels  to  the 
levels  measured  at  the  entry  port  of  the  waveguide.  It  is 
considered  that  the  following  discussion  supplies  that 
relation. 

c.   Total   Power  at  Waveguide  End 

The  expression  for  P  (x)  in  (V.  C.  3)  above 
gives  tae  spatial  distribution  of  power  as  a  time  average. 
The  integral  of  that  expression  over  the  face  of  the 
waveguide  gives  the  total  time  average  power  in  watts 
delivered  to  the  open  end  of  the  waveguide  (P  (o) ) : 
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P"(o)  =  P        =  f   /IP  (x)ldxdy 
total     o  o 
avg 

a 
=  f   K  sin2(Trx\dx       K  =  jStMLa^lHoJ^b 

P(o)  =  4^ua2|.Hol£b  (V.  C.  4) 

If  the  power  measurement  system  were  lossless  and  there  were 
no  reflected  power,  the  power  measured  as  going  down  the 
waveguide  (P  (s) )  would  be  the  value  F(o).  By  comparison 
with  P  in  (V.  C.  3) ,  it  is  seen  that  the  power  being 
delivered  at  the  center  of  the  waveguide  (P  peak)  is: 

P(x)  peak  =  2   P  (o)  £  (V.  C.  5) 

a73 

The  relationship  now  to  be  investigated  is  that  between  the 
total  power  going  down  the  waveguide,  measured  at  the  power 
measurement  system,  V (s) ,  to  the  total  power  being  delivered 
at  the  open  end  of  the  waveguide,  F(o). 

d.   Relation   between   Total   Average    Power   at 
Waveguide  Source  and  Load 

Because  of  reflections  from  the  load  there  are 
standing  wave  patterns  in  the  waveguide.  It  is  known  that 
the  current  and  voltage  standing  wave  patterns  are  periodic 
with  a  period  of  half  the  waveguide  wavelengths  and  are 
shifted  with  respect  to  each  other  by  a  quarter  wavelength. 
Figure  15  shows  the  patterns  for  an  open  ended  transmission 
line. 33  Since  the  power  delivered  is  a  product  of  voltage 
and  current,  it  would  be  expected  that  there  be  regions  of 
maximum  and  minimum  power  availability.  In  a  real  system, 
if  there  is  no  load  present  to  accept  that  power,  then 
currents  shall  flow  on  the  outer  surface  of  the  waveguide, 
and  the  power  shall  be  expended  in  the  losses  of  the 
waveguide. 

"Ibid^,  pp.  219  -  220. 
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Figure  15 
The  presence  of  these  power   standing   waves  can   be   shown 
analytically: 

Defining  the  voltage  reflection  coefficient   to 
be  the  ratio  of  reflected  to  forward  electric  fields 

and  the  origin  of  the  z-axis  to  be  at  the  film   end   of   the 
waveguide  as  shown  in  Figure  16: 


-1 


Rl 


Figure  16 
Wave  guide  Transmission  Line  Coordinates 

At  a  point  z  =  -£  the  incident   electric   field 


VF  is 


V_  =  V  e  =Ver 

^     o  o 

where  V   is   the   peak   value  (in   time   and   space)  of  the 

o  _^ 

incident  electric  field.   The  reflected  field  VR  is 


VR  Mvje 


•   =  P  V  e   r 
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and  the  total  electric  field  at  point  z  =  -J!  is 

Vp  +  V^  =V  (e  r  +  f  e  r  ) 

or  generally 

~j£z     +  j£z 

V  (z)  =  V  (e   '  +  f  e  r    ) 

-j/Sz        H    jftz 

I(z)    =   Vo    (e      r  -  fe    r  ) 
Zo 


(V.  C.  6) 


the  time  average  power  density  at  any  point  z  then  is 
P(zJ  =  J.  Fe  {V  (z)  »I  (z)  } 

-iH.fTo.[i  *re^z  -rYj*"    -rrV*} 

2   (E5~ L  J    J    (V.  c.  7) 

where  we  indicate  the  complex  conjugate  of  a  function  A  by 
A   and  take 


V(z,t)  =  V(z)  cosut 

I(Z,t)   =  I(Z)   COS  (6jt  +  0 


) 


V    is   a   real   number,  the  peak  value  of  incident  electric 
o 

field,    0  is  a   time-phase  difference  between  electric  and 
magnetic  field  peaks. 

At  the  source,  where  the  power  measurement 
setup  is  located,  the  power  measured  must  be  V (z)  evaluated 
at  z  =  -  2  ,  which  must  also  be  the  power  being  measured  at 
the  source  after  corrections  are  made  for  losses  in  the 
measurement  system.  That  is,  P (s)  is  the  resultant  of 
incident   and  reflected  power 

V(s)    =   *F-  P* 

=  J.  Ee  f  Vo£  n  +  (Te-J2fs    -p*eJ2fl)  -  pri**|  e-j*l 

He  can  relate  P(s)  to  P(o)  by  noting  that  by   (V.  C.  7),  for 
z=0 


(o)  =  I  He  f  Vo£  (1+  ^-^*^pp*)  e'J  7 


so  that 
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m =  f(p) 

P(o)  =  f  (p)P(s) 

P(o)  =  Ref n  ,.{  r.7  r  *)  .-  r,      7  .  p<S)     (v.  c.  8) 


>f   1  +  /  r-  r  *)  -  rr*         i   .  p<: 


e.  Relation  between   Total   Power   at  Source   and 
Power  Measured  by  Meters 

The  final  step  in  finding  the  value  of  power 
being  delivered  is  to  relate  P (s)  to  the  power  levels  read 
at  the  meters.  The  losses  in  the  power  meters  and 
thermistor  mounts  are  provided  by  the  manufacturer  as  a 
parameter  called  the  calibration  factor  (CF)  ,  which  compares 
meter  reading  (PM)  to  the  power  incident  on  the  thermistor 
mount  (Px  ) 

CF  =  P* 
P 

For  the  general  case  in  which  there  is  both  forward  (PF)  and 
reflected  (PR)  power  being  read, 

P  =  PF  -  PR 

=  IW_  -  Pm%_ 
CFF    CFR 

where  PKf  ,  PMR  are  the  forward  and  reflected  power  meter 
reading  and  CFf ,  CFR  are  their  respective  calibration 
factors. 

In  these  experiments  the  power  levels  were  such 
that  in-line  attenuators  RF  and  \    were  required.   Thus, 

pcs>  =    [ir_  ♦  Rf]     -    [^  ♦  «,]  (V.  C.  9, 

f.  Summary 

In  summary,  beginning  with   assumptions  of  TE 
mode   operation   above   cutoff   the   spatial  distribution  of 
average  power  across  the  waveguide  (P  (x) )   was   found.   Next 
the  spatial  total   average   power   at   the   load   (F(o))  was 
developed.   That  value  was  then  related  to  the  average  power 
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applied   at  the   source  (P"(s)).    Finally,  the  relationship 

between  the  source  power  P(s)  and  the  incident  and  reflected 

power  meter  readings  (PM(t,PMR)  was  established.     Thus  the 

equations   used   to   relate   meter   readings  to  peak  spatial 

_*. 
average  power  P{x)  evaluated  at  x  =  a  were: 

p"(x)peak    =2     P  (o)  £  (V.    C.    5) 

aE 

F(o)    =    Re        1    ±__r_r._V*  .-  fT  *  .    F  (s)  (V.    C.    8) 

1  ♦  rfm   -rV*Pr  -rr* 

P(s)    =    Pr         -   Pf 

=  [&;-  *  Efl  -  Clr- *  E*]  (V-  c-  9) 

where 

a  =  x-dimension  of  waveguide 

P  =  voltage  reflection  coefficient 

Ji.  =  length   of   waveguide   from   lead  to  power  measurement 
setup 

PH  =  power  meter  reading 

CF  =  thermistor  mount  calibration  factor 

R  =  in-line  calibrated  attenuators 
These  formulas  were  used  to  measure  the  power   delivered   to 
the  center  of  the  end  of  the  waveguide. 

2 .  Determination   of Power. ^Delivered hl^& Reflector 

System 

A  driven  monopole  and  60°  corner  reflector  over  a 
ground  plane  were  used  to  produce  the  field  for  some  of  the 
experiments.  The  use  of  a  conducting  reflector  permitted 
application  of  image  and  array  theory  to  determine  the  field 
in  the  vicinity  of  the  system.  Such  a  determination  was 
considered  necessary  in  order  to  select  the  placement  of  the 
reflector  and  of  the  target  element  which  would  optimize  the 
field  at  the  target  element.  The  geometry  of  the  system  is 
shown  in  Figure  17. 

The  field  reflected  by  a  perfectly  conducting 
ground  plane  and  corner  reflector  can  be  considered  to  be 
produced   by   five   image   elements   arranged  in  a  hexagonal 
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pattern.34  Phasor  addition  of  the  individual  fields  from 
all  six  elements  then  gives  an  approximation  of  the 
resultant  field.  Because  of  the  different  element-to- target 
distances  involved,  there  are  phase  and  magnitude 
differences  which  may  vary  significantly  along  the  axis  of 
symmetry.  The  following  development  explores  these 
possibilities. 


£>0°  corner 
Elector 


Qheun 


Figure    17a 
Perspective   View   of   Reflector    System 


Key: 

A    Driven  eUvnent 

n    Trv»aeie   Clement 
O  Tdr-aet  clement 


a* 


Figure  17b 
Dimensions  of  Reflector-Image  System 


34Kraus,  p.  330 
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a.   Fields  of  a  Single  Dipole 

The  near  field  of  a  single  dipole  has  been 
derived35  giving  the  following  expressions  for  the  electric 
and  magnetic  field  components  of  a  quarterwave  monopole  over 
a  conductive  ground  plane.  The  subscript  m  shall  be  used  to 
indicate  that  the  rath  element  is  being  considered. 


/   N 


-DPR.*      -j£**m     „       -J^* 


E*   =  Ec   z  =  -jEfc  '  e_   ♦  e_      -  2cosftHe_ 


(V.  C.  10) 


m  0  •  t 


m 


E*0'  =  30lr»' 

where  Im'  is  the  element  source  current,  H  is  the  monopole 
length  (dipole  half-length)  ,  &  is  the  propagation  constant, 
n  is  377  ohms,  z  and  96  are  unit  vectors  in  their 
respective  directions,  and  the  Rtm  ,  rm,  and  ym  are  defined 
for  the  mth  element  of  Figure  17b  by  the  geometry  of  Figure 
18.   From  that  figure  it  is  seen  that 


rm=V*m  +  z2 


These   variables  can  be  reduced  to  multiples  of  a  wavelength 
(X)  by  defining 


H  =  hX 
z  =  s\ 


d  =  C* 


35Jordan  and  Balmain,  p.  336. 
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Figure  18 
Geometry  of  a  Monopole  over  a  Perfectly  Conducting 

Ground  Plane 

Substituting  (V.  C.  11)  in  (V.  C.  10)  results  in 


E    '    =    -jE 


*Q 


fcos8im  +  ,cosGam    -    2cos8s  cosQ<tm? 

[I       rim  r*m                         ?m                J 

-  ifsinSim  +    sir.8z.Tw     -    2cos&3  sin&sm ) 

L       rim  rtm                       f>m                 J 


where 


On*  =  2Trr1M 

02m=  2>r6m 

Q*    =  2irh 

G*™  =  21ffm 

E^       «  E»  ' 


Separating    E    •    into    its   real    and    imaginary    parts  yields 

m 


Ee{E^'}    =   -E^      sin©2_rn    +    siH^S.  ~    2cos0asin6gm 


0        rn 


im 


P« 


(V.    C.    12a) 


Im{E,J}      =    -E,        COs(W      +     COS©Rrn     -     2cOsD3_COS04m 


o   r, 


it* 


?* 


A  similar  operation  on  H    results  in 

m 

Re£H0»vU     =  _£_Sfi_  (sinGim    +    sin©tm  -    2cos©s  sin0nm  ) 


'm 


(V.    C.     12b) 


Im(*W  )      =     Ey        (COsOim     +     COsOim     -     2cOS©3  COS^Hth  ) 
7  Km 
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From  the  geometry  of  Figure  17b  it  is  clear  that  all  the  lit 
shall  not  be  aligned,  so  that  the  magnetic  field  resulting 
from  all  six  elements  shall  have  two  components,  one 
parallel  (Hy)  and  one  normal  (Hx)  to  the  axis  of  symmetry. 
These  components  are  determined  by  the  geometry  of 
Figure  19. 


A 


Figure  19 
Normal  and  Parallel  Components  of  the  mth  Magnetic  Field 

From  that  figure  and  the  law  of  cosines,  the  included   angle 
(<*)  is 

oc  =  cos-*  Yx  2  +  ym2  +  2dy,t, 
2  (3  +  y,  TyZ 

=  COS'i  K^2_+_Km£_+_2CK1_ 
2(C  +  KjJK^ 

so  that  HtfnJ  can  be  broken  into  its  components, 

*W    =   -    CHxm«    x   ♦    H-    y] 


H^    = 


H0rn«     cos«( 


Hynl    =    H^    sin* 

Examination  of  the  geometry  of   Figures   17b   and   19   shall 
reveal  that  the  Hy  components  of  elements  2  and  3  cancel  the 
Hy  components  of  elements  6  and  5  respectively, 
b.   Effects  of  Reflection 

The  development  to  this  point  assumes  a 
perfectly  conducting  ground  plane  and  corner  reflector  with 
perfect  connection  to  the  ground  plane.  In  such  a   case   the 
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reflection  coefficient  is  -1.  The  (tangential)  electric 
field  is  reversed  upon  reflection.  The  tangential  magnetic 
field  component  is  unchanged,  while  the  normal  magnetic 
field  component  changes  sign.  The  net  effect  is  that  of  the 
source  current  in  the  image  element  being  180  degrees  out  of 
phase  with  that  of  the  driven  element.36  Ifr  however,  the 
reflector  is  not  in  perfect  contact  with  the  ground  plane 
the  reflection  coefficient  (T)  becomes  complex.  In  such  a 
case  the  image  source  current  suffers  a  phase  lag  as  well  as 
attenuation,  and  the  reflected  electric  field  (E.^)  is  the 
product  of  the  reflection  coefficient  and  the  incident  field 
(E^m*)  for  each  reflection: 

P  =  iPle 


m    '    m 


For  multiple  reflections  the  resultant  electric  field 
arriving  at  the  target  element  is  the  original  source  field 
multiplied  by  the  appropriate  power  of  the  reflection 
coefficient 

fU  =  1  for  m  =  1  (V.  C.  13a) 

ft„  =  P  for  m  =  2,4,6 

T*  =  (rj  2  for  m  =  3,5 

and 

C*  =  r«£*  =  ru^  £  +  Hym«  y) 

and  (V.  C.  13b) 

c.   Summation  of  Component  Fields 

For  the  array  of  real  and  imaged  dipoles  shown 
in  Figure  17b  at  a  point  yt  =  K£X  along  the  axis  of  symmetry 
from  the  driven  element  the  total  fields  are  found  by  phasor 
addition  of  the  contributions  from  each  element,  taking  into 


36Jordan  and  Balmain,  p.  470  -  471 
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consideration   the   phase   and  angular   variations   of   the 
electric  and  magnetic  fields  described  above: 


where 


E*t  =  E*t  z 


E^    =    Re{E^}    +    jIm{Ej    } 


(V.    C.    14a) 


Re{E?t}    =  E    (Re{E,,J) 


m=1 


Also, 
where 


xt 


Im{E«}     =ra?1(Im{E*^) 


H0t    -    ~    HKt  x 


(V.    C.     14b) 


m=  i  m  = 1  m=  1 


*hv 


and   the    mth   contributions   to    these      sums      are      described      by 

equations    (V.      C.     12)    and    (V.      C.     13). 

d.       Total   Time  Averaged    Power    Density 

The        total      time      averaged      power      density      ?t 

impinging   on    that    point    y         is      then         found        by        applying 

i 

equation  (V.  C.  2)  to  these  total  fields: 

P    =  ±   Re  {%  } 
av    2 


_».  t,     -*. 


where 

p*  =  E«x  *W 
That  cross  product  is  found  from 

x    y    z 


(E^x  H^c)  =  det 


Ex    Ey   Ez 
-HX*   0     0 


=  ft  =  p,t  y  ♦  pw  i 


so  that 


P^t  =  -EzHx 
P^  =  EyHx* 


(VI.  C.  15) 
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From  the  geometry  of  Figure  18  R^m  can  be  found  in 
terms  of  z  and  the  various  y^,  hence  r^  in  terms  of  s  and 
Km: 


r   2  =  f^i.tn\2    =     (Km*  +  s2  +  h2)  -  2hs 

(V.  C.  16) 
r   2  =  (Rh±\ 2  =  (Km2  +  S2  +  a*)  +  2hs 

2m      ^     \    J 

To  find  the  Km  we  use  the  geometry  of  Figure  17b  and  the  law 
of  cosines,  where  the  angles  oc  and  |3  are  seen  to  be  sixty 
and  120  degrees  respectively: 

y2=(d  +  y)2+d2-  2  (d  +  y  ) •d»cos60° 

2  1  1 

y2=(d    +    y)2    +    d2-    2(d   +    y   )  «d»cos120° 
When    simplified    these   give 


y2=d2+yd+y2 

2  11 

y    2    =    3<}2    +    3y    d    +    y    2 

3  1  1 


(V.     C.     17) 


Taking 


results    in 


By    symmetry,    then 


y,  =  \x 

y      =   K    \ 

2  2  (V.      C.       18) 

y      =   K    x 

3  3 

y      =  K    x 
*  ♦ 

y     =  K   X 

5  5 

y      =  K    X 

6  6 

d    =   C  X 


K    2    =    c2    +    CK      +    K    2 

2  11 

K    2    =    3C2    ♦    3CK      +    K    2  (V.    C.     19) 

3  11 


K      =    2C    +    K 

4  1 

K      =    K  (V.     C.    20) 

5  3 

K       =    K 

6  2 


Then,       using      equations       (V.    C.    19)         and       (V.    C.    20)         with 
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equations  (V.  C.  14)   the  total  power   present   at   a   point 
y   =  K  X  along  the  axis  of  symmetry  can  be  found  by  summing 
the  contributions  from  all  six  elements, 
e.   Measurement  of  Power 

The  instrumentation  of  Figure  12  was  used  to 
monitor  the  power  being  fed  to  the  driven  element  of  the 
system.  Additionally  a  probe  was  used  to  measure  the  field 
along  the  axis  of  symmetry.  The  probe,  shown  in  Figure  28, 
sampled  the  radial  component  of  power  flow.  The  sample  was 
then  measured  with  a  thermistor  and  power  meter  arrangement 
such  as  that  described  on  page  37  above.  Since  the  probe  is 
aligned  vertically,  it  is  considered  that  it  only  sampled 
the  radial  component  of  power  flow.  In  order  to  compare 
experimental  results  with  the  expectations  of  equation  (V. 
C.  15)  ,  the  resultant  value  of  P\jt  was  first  integrated 
numerically  using  Simpson's  rule  over  the  range  in 
wavelengths,  which  corresponded  to  the  length  of  the  probe: 


E*t 


.05X         r.05X         f.05X 
=  j     Pr  ds  =  J     Px  ds  +  j     Py  d: 
0  0  0 


A  listing  of  the  program  written  to  carry  out 
the  indicated  summations,  products,  and  numerical 
integration  is  contained  in  Appendix  C.  It  was  found  that 
fluctuations  occured  in  the  power  being  delivered  to  the 
system  (P(wc  )  when  either  the  probe  or  corner  reflector  was 
in  the  very  near  field  of  the  driven  element.  The  effect  of 
these  variations  on  the  measured  power  magnitude  (P4tm  )  was 
corrected  for  by  applying  a  correction  multiplier  ( §  )  such 
that 

The  incident  power  was  measured  with  the  offending  element 
moved  far  away  from  the  driven  element.  This  value  (P0 )  was 
taken  as  a  reference.  The  ratio  between  distant  spacing  to 
near  field  power  values  was  taken  as  the  multiplier,  that  is 
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4-1 
K 


Figure  22 
Power  vs  Spacing  in  the  Field  of  a 
60°  Corner  Reflector  Radiator 
(Monopole  to  Probe  Distance  is  Principal  Axis) 


Figure  23 
Power  vs  Spacing  in  the  Field  of  a 
60°  Corner  Reflector  Radiator 
(Corner  to  Monopole  Distance  is  Principal  Axis) 
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Figure  24 
Power  vs  Spacing  in  the  Field  of  a 
60°  Corner  Reflector  Radiator 
(Monopole  to  Probe  Distance  is  Principal  Axis) 


Figure  25 
Power  vs  Spacing  in  the  Field  of  a 
60°  Corner  Reflector  Radiator 
(Corner  to  Monopole  Distance  is  Principal  Axis) 
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Figure  26 
Normalized  Power  vs  Spacing 


Figure  27 
Normalized  Power  vs  Spacing 
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I  ■  «t- 


T*C 


Figures  20  and  21  show  plots  of  the  results  of  integration 
for  different  values  of  the  parameter  T  .  Figure  22  through 
25  show  the  composite  dependence  on  spacing.  Figures  22  and 
23  give  power  on  a  dbm  scale.  Figure  24  and  25  show  power 
in  milliwatts.  Figures  26  and  27  are  normalized  to  show  the 
shift  in  peak  power  with  different  spacings. 


Figure  28 
Field  Power  Sensing  Probe 

3 •   Determination  of  the._Ef  f  ects.. of  System  Resonance 
a.   Deter mimation  of  Film  Impedance 

To  determine  the  effects  of  system  resonance  it 
was  considered  necessary  to  model  the  film,  then  determine 
the  response  of  that  model  to  an  applied  field.  A  resonant 
circuit  was  considered  likely; 


^ 


£ 
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This      is    seen    to   be    equivalent: 


"C 


T 


The  capacitance  is  provided  by  the  gap  between  waveguide  and 
the  conductive  film  negative.  The  resistance  and  inductance 
are  present  in  parallel  in  the  conductive  media  (the  film 
negative  and  the  activating  chemicals) . 

The  Laplace  transform  of  the  impedance  of  such 
a  circuit: 

z  (s)  =  _L_  + * 

SC    -i  +  i 

R   Ls 

=  S2RLC  +  SL  +  R 
S^LC  +  SRC 


+  /L  -  Rzcl.  1  S  !  L  -  P2C 
[   LC    J   [s2  +  s(k/l)l 


Z  (s)  =  R  +  f  L  -  R2C 


letting  a  =  P,  b  =  L  -  R?C,  c  = 

LC 

becomes 


,  d  =  JL  the  impedance 
L  -  R2C       L 


R 


Z(s)  =  a  +  b(s  +  c) 
s(s  +  d) 


(V.  C.  21) 


b.   Determination  of  System  Time  Response 

The  experiment  consisted  of  applying  a  signal 
to  the  film  to  determine  its  characteristics,  and,  later,  to 
obtain  images.  The  response  of  the  model  to  such  a  signal, 
considered  to  be  cosinusoidal  and  of  unit  amplitude  is  found 
from  a  transmission  line  analogy,  using  Laplace  transforms. 
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i (t)  =  cosot 
I(s)  = s 

S2  +  UJ2 

V(s)  =1  (s)  •  Z(s)         (V.  C.  22) 

By  expanding  (V.  C.  22)  using  partial  fractions  and 
combining  certain  terms,  the  total  response  transform 
becomes 

V(s)    =    sLa^^Cj^lcos^cu^JJ    +   _LbCi._M  sinjjfcupVlJ  -    bCaJ.^) 
s2    +  u>2  s2    +  ^"2  s*  +   "5 

where 


c 

(w) 

1   k? 

+   02 

1 

wyd2 

+  CO  2 

c 

/-c-i- 

d     ■ 

d2      +    w2 

0(a>)    =   tan-i  {£}    -   tan-i  {£} 

for    which   the   inverse  transform    is 

-dt 
V(t)     =    K     y  cos[wt    -    9(w)  ]   -    K     (o>)  e  (V.    C.    23) 

1  2 

WhGTG  — — — — — — — — —  

K    (u>)    =V  a2    +    b2c    z  (u>)     +    2abC    (u>)  cos[0  (o)  ] 
K    (io)    =    bC    (co) 

2  2 


e<«)  = 


tan~1  fJ'iuWV^MJ  n) 

(  a  +  DC1  (oj)  COS|_0(cO)  ]  j 


D.    DETERMINATION  OF  THERMAL  PARAMETERS 

Various  thermal  parameters  were  measured.  All  thermal 
measurements  were  made  with  a  thermocouple  and  a  digital 
thermocouple  indicator37.  The  ambient  temperature  was 
22°  Centigrade  (71.6°F).  Variations  of  temperature  with 
time  were  determined  under  various  conditions,  and  changes 
in  film  density  were  determined  for  different  initial 
temperatures. 

370MEGA  Thermocouple  Indicator,  Model  DS500,  Type  J 
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1-       Cqoling->Box_Te^perature 

a.  Cooling  Rate  of  Box 

It  was  considered  necessary  to  determine  the 
time  required  to  reach  a  steady  temperature  inside  the 
cooling  box  (Figure  5)  in  order  to  ensure  uniform  data.  The 
box  was  filled  with  the  appropriate  coolant.  The 
thermocouple  was  passed  through  a  small  hole  in  the  side  of 
the  cooling  box  at  the  level  at  which  the  film  would  be 
placed.  The  top  to  the  box  was  put  in  place  at  time  zero 
and  the  temperature  was  then  recorded  at  5  second  intervals 
until  steady  state  had  been  reached.  This  experiment  was 
repeated  several  times,  and  the  average  temperature  for  each 
time  point  obtained  and  plotted  versus  time. 

b.  Effect  of  Film 

The  temperature  changes  inside  the  box  when  a 
f-ilin  packet  was  inserted  was  similarly  measured  and  plotted. 

2.  Ac tivat or _ Chemical  Freezing-point 

An  experiment  was  conducted  to  find  the 
approximate  freezing  point  of  the  chemical.  A  thermocouple 
was  frozen  into  a  portion  of  the  chemical  and  readings  were 
taken  as  the  chemical  warmed  to  ambient  temperature.  The 
freezing  point  was  considered  to  be  the  temperature  at  which 
the  warming  curve  was  minimum,  since  at  that  point  energy 
was  being  used  to  change  state  from  solid  to  liguid. 

3 .  Temperature Variations in Film Due   to Cooling 

Method 

Two  methods  were  used  to  pre-chill  the  film. 
Experiments  were  conducted  to  determine  the  temperature 
changes  in  the  film,  with  time  for  the  different  cooling 
methods,  and  thus  provide  a  basis  for  determining  the 
initial  temperature  of  the  film  when  first  placed  in  the 
field. 

a.   Dry  Ice 

Dry  ice  was  used  to  cool  the  film  to 
temperatures  between  -30°  and  +20°C.  To  ensure  consistent 
results  the  cooling   box   was   always   filled   to   the   film 
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opening  with  dry  ice.  Crushed  dry  ice  was  used  to  fill  all 
niches.  The  thermocouple  was  inserted  through  the  box  wall 
and  used  as  a  film  support  for  all  experiments  where  the  ice 
was  used.  The  probe's  position  resulted  in  a  fairly 
constant  cooling  position  for  the  film.  Also  the  initial 
temperature  of  the  box  before  film  insertion  could  be 
monitored.  Consequently  the  experiment  could  be  conducted 
under  constant  conditions.  The  temperature  at  the  probe 
consistently  was  50°C  before  film  insertion, 
b.   Liquid  Nitrogen 

For   cooler   temperatures   liquid   nitrogen  was 
used  as  the  coolant.   N   was  poured  into  a  pre-cooled  box  30 

2 

seconds  before  film  insertion.  Precooling  was  accomplished 
by  pouring  60  ml.  of  the  coolant  into  the  box  1  minute 
before  film  insertion  and  removing  it  45  seconds  before  film 
insertion. 

To  measure  the  temperature  inside  the  film  this 
procedure  was  modified  slightly.  The  film  was  inserted  in 
the  cooling  box  at  time  zero.  The  temperatures  were  obtained 
with  the  thermocouple  previously  inserted  through  a  small 
(2mm  square)  hole  into  the  space  between  print  and  negative 
in  the  film  packet.  The  temperatures  at  5  second  intervals 
were  recorded  and  averaged  over  several  runs  for  each  point 
in  time,  then  plotted  versus  time. 
4 .   Warming  of  Film  on, Removal 

Finally,  data  was  taken  to  determine  the  rate  at 
which  the  film  returned  to  ambient  temperature  after  removal 
from  the  cooling  box. 
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V I .   PfiESE NT AT 10 N_  A  N D_E V  A LUA T 10 N _0 F_E X P E RIME N T A L_R E S UL T S 

A.    SUHBARY  OF  PARAMETERS 

In  order  to  properly  consider  the  results  obtained  it 
is  considered  useful  to  summarize  the  various  process 
parameters  and  their  expected  relationships. 

1 .   Principal...  parameters 

The  process  may  be  viewed  as  a  system  for  which 
the  input  is  applied  energy  and  the  output  is  image  density. 

a.  Applied  Energy 

The   energy   is   applied  to  the  system  in  three 
stages: 

(1)  In  the  sensitization  preexposure  phase, 
the  parameters  are  the  amount  of  preexposure  to  blue  and 
green  light. 

(2)  In  cooling  the  film,  the  principal  energy 
parameter  is  the  film  packet  temperature,  which  is  itself 
dependent  on  the  temperature  of  the  pre-chill  box. 

(3)  In  the  radiation  phase,  the  critical 
factors  are  initial  temperature  of  the  film,  ambient 
temperature  in  the  room  where  the  experiment  is  being 
conducted,  thermal  response  of  the  chemicals  involved, 
strength  of  the  applied  energy  field  and  duration  of 
radiation. 

b.  Image  Density 

The  density   (D)   of   the   resultant   image   is 
dependent  on  the  interaction  of  all  the  above  parameters. 
2-   Experimental _0bjec five 

The   overall  experimental  objective  was  to 

a.  Characterize  the  various  parameters. 

b.  Characterize  the  image  density  response  to 
those  parameters. 

c.  Examine  that  response  in  order  to  determine  an 
optimal  set  cf  parameters  which  would  provide  the  best  image 
of  an  electromagnetic  field. 
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These  objectives  were  to  be  met  by  three  sets  of 
experiments.  The  first  set  may  be  considered  to  be  a  set  of 
static  tests,  where  a  single  parameter  is  varied  on  a  given 
run.  The  results  of  these  tests  are  presented  in 
paragraphs  B  and  C  below. 

On  the  basis  of  these  results  certain  preferred 
sets  of  parameters  were  selected  with  which  the  waveguide 
experiments  were  to  be  conducted.  This  set  of  experiments 
may  be  considered  to  be  a  set  of  dynamic  tests,  since  these 
measure  the  dynamic  qualities  of  LB  and  A.E,  that  is,  what 
variations  in  film  density  (i.e.  contrast,  AD)  occurs  in 
response  to  a  spatial  variation  in  field  strength  (AE)  . 
These  experiments  are  discussed  in  paragraph  C. 

The  final  set  of  experiments  is  the  confirmation 
set,  which  were  conducted  to  confirm  the  conclusions  drawn 
on  the  first  two  sets. 

B.    THERMAL  VARIATIONS 

1 •   Cooling  Box_Temperature 

With  the  cooling  box  filled  with  crushed  dry  ice 
to  the  film  entry  port,  the  top  was  put  in  place  at  time 
t  =  0.  The  measurements  were  taken  at  the  position  of  the 
film,  as  seen  in  Figure  29. 


Cooling  Boy  n. 


TWevrmo  coo  pie  Lead 


_":  -._-       £  ( 


^r«fcV»c<it?irtkXce. 


Figure  29 
Measurement  of  Cooling  Box  Temperatures 
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The   results  shown  in  Figure  30  show  that  the  cooling  box  is 

effectively  at  steady  state   4  minutes   after  being   loaded 

with   dry  ice.  The  steady  state  reached  after  15  minutes  was 
consistently  within  a  few  degrees  of  -50°C. 
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Figure  30 
Temperature  versus  Time  for  Cooling  of  Chill  Eox 

Figure    31   shows   the   change   which   occurs   in   the   box 

temperature  when  a  film  packet  is  inserted  at  time  t  =  0. 
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Figure  31 
Temperature  versus  Time  for  Warming  of  Chill  Box 
When  Film  Is  Inserted 
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2-   Film  Temperature  Variations 

Figure  32   shows  the  method  used  to  record   changes 
in  the  film  temperature  during  the  prechill  process. 
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Figure  32 
Measurement  of  Film  Temperature  Changes 

Figures  33  and  34  show  the  curves  for  dry  ice  and 
liquid  nitrogen  cooling,  respectively.  These  data  are  the 
basis  for  initial  film  temperature  values  assumed  throughout 
for  other  experiments  and  data.  These  curves  are  of  the 
form 


-t/y  "t/  y 

Ji  +  K  e    J 


T(t)  =  K   +  K  e 


TT  has  been  seen  to  generally  be  of  the  order  of  50  seconds, 
l 

with   ~r  approximately  180   seconds.    K  ,   K  ,  and   K   have 
J2   r  r  J  12  3 

typical  values  of  -35,  7 ,  and  50  respectively  for  dry  ice, 
75,  70,  and  27  for  liquid  nitrogen.  Each  datum  represents 
an  average  over  3  runs  (CQK)  and  5  runs  (Ni).38  The  average 
standard  deviation  for  these  data  points  are  2°. 4  and  2°.  2 
for  the  dry  ice  and  liquid  nitrogen  curves,  respectively. 


used. 


38See  Appendix  D  for  definitions  of  statistical  methods 
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Figure  33 
Film  Cooling  with  Dry  Ice  Coolant 
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Figure  34 
Film  Cooling  with  Liquid  Nitrogen  Coolant 
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The  film  reached  apparent  steady  state  after 
approximately  five  minutes.  Approximately  five  minutes 
later  the  film  was  removed  from  the  cooling  box  and  placed 
upright  in  a  position  similar  to  that  being  used  to 
irradiate  the  film.  Figure  35  shows  the  time  response  of 
the  film  to  these  conditions.  Time  zero  represents  the  time 
at  which  the  film  was  removed  from  the  box.  This  curve  is 
of  the  form 


+  K  e    Ji  +  K  e    « 


T(t)  =  K 

12  3 

fiepresentative      numbers      for      K    ,    K    ,     -v      ,    K    ,    and      y     are, 

12      13-2 

respectively  22,5,15  seconds,  55  and  95  seconds. 
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Figure  35 
Temperature  versus  Time  for  Film 
when  removed  from  Cooling  Box 

The   data   of  Figure  35  was  averaged  over  4  runs,  and  had  an 
averaae  standard  deviation  of  1°.4.39 


3 9  ibid. 
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Figure  36 
Warming  Curve  of  Activator  Chemical 

* ^ .   Activator  Chemical  Freezing  Point 

Figure  36  shows  a  one  run  result  of  an  experiment 
to  find  the  approximate  freezing  point  of  the  chemical.  The 
minimum  slope  (D2T  =  0)  occurs  when  the  temperature  reaches 
-7°C#  which  is  taken  as  the  freezing  point  of  the  chemical. 

C.    DENSITY  VARIATIONS 

1 .   Density_ Variation  .with  Sensitizing  Pre-exposure 

Figure  37  and  38  display  the  density  variations  as 
a  function  of  two  parameters  on  logarithmic  axes.  These  are 
the  results  of  the  experiment  described  in  paragraph  V.  A. 
above.  The  horizontal  axes  represent  the  composite  amount 
of  pre-exposure  to  filtered  blue  and  green  light  (Log  H 
(Blue)  and  Log  H  (Green)).  The  vertical  axis  is 
representative    of   the   average   film   density   (D)   which 
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resulted  from  such  pre-exposure.  Each  figure  provides  a 
view  of  the  same  data  from  different  viewing  angles  to  allow 
observation  of  different  aspects  of  the  film's  response 
characteristic.  These  figures  reveal  a  significant 
variation  in  dynamic  range  over  the  field  of  measured 
exposure  levels.  The  steepest  slope  appears  to  exist  in  the 
vicinity  of  (-1.55,  -.5)  which  corresponds  to  a  blue^green 
pre-exposure   of   approximately    3  sec,    blue.  _lsec.  green. 

r  100  10 

The  gradient  appears  to  be  aligned  diagonally  to  both 
horizontal  axes.  The  nearest  corner  of  the  density  function 
surface  shows  the  lowest  response,  corresponding  to 
preexposures  of  3  seconds  for  blue  and  2  seconds  for  green 
light. 


Time  (%tc} 


Figure    37 
Average   Film   Densitv    versus    Preexposure   Time 
(Principal    Axis' is    Blue   Preexposure) 
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Figure  38 
Average  Film  Density  versus  Preexposure  Time 
(Principal  Axis  is  Green  Preexposure) 

2 •   Density,  Variation^ with  Initial_Temperature 

Figure   39  and  40  show  the  changes  in  density  which 

occur   for   different   values   of   initial   temperature   and 

applied  energy. 

The  steepest  slope  appears  to  be  in  the  vicinity  of 

the  temperature-energy  coordinates  of  (0°,  1000  joules) . 

The  gradient  is  aligned  diagonally   to   both   axes, 

showing  an  overall  dependence  on  total  energy  applied  during 

the  radiation  phase. 
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Figure  39 

Density  versus  Radiation  farazeters 

(Principal" Axis  is  Nominal  Radiated  Energy) 
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Fiaure  40 

Density  versus  Radiation  Parameters 

(Principal  Axis  is  Initial  Temperature) 
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D.  DYNAMIC  RANGE  MEASUREMENTS 

The  waveguide  experiments  were  conducted  to  determine 
the  dynamic  range  of  the  film,  to  be  characterized  as  the 
product  &  of  contrast  (AD)  and  change  in  energy  (AE) .  These 
experiments  were  hampered  by  a  factor  not  previously  noted: 
thermal  noise,  that  is,  unwanted  images  due  to  ambient 
temperature  in  the  room.  The  static  tests  were  conducted 
with  a  uniform  image  (AD  =  0)  so  the  incremental 
contribution  of  ambient  temperature  was  not  seen.  However, 
the  waveguide  experiment  required  simultaneous  imaging  of 
high  and  low  energy  regions.  An  expected  variation  of 
density  occured  in  response  thereto.  Additionally,  though, 
there  was  a  uniform  density  distribution  due  to  the  uniform 
ambient  temperatures  near  the  experiment.  The  effect  was 
minimized  by  using  long  pre-exposure  times,  e.g.  2  seconds 
blue,  1  second  green  and  low  initial  temperatures,  e.g. 
-30°.  This  however  restricted  operations  to  the  flatter 
regions  of  the  surfaces  presented  in  Figures  33  and  35. 
These  flatter  regions  have  a  much  lower  and  generally 
uniform  dynamic  range  as  measured  by  (T   =  adae. 

The  phenomenon  is  directly  analogous  to  that  of  a 
receiver  or  control  system  in  which  high  sensitivity  and 
dynamic  range  is  accompanied  by  poor  performance  because  of 
noise  response,  so  that  optimal  performance  is  ultimately 
found  to  be  derived  from  a  raised  a.g.c.  threshold  or  a 
damped  system.  Figures  41  and  42  show  two  photographs,  one 
taken  at  a  high  sensitivity  point  on  the  surfaces,  the  other 
taken  at  a  low  sensitivity  point.  The  pref erability  of  the 
latter  is  obvious. 

E.  FIELD  MEASUREMENTS 

Figure  43  shows  the  normalized  precoraputed  and 
normalized  experimental  values  of  relative  power  in  the 
field  of  the  dipole  and  reflector  system,  as  a  function  of 
the   corner^to-dri ven-element   distance   for   a    reflection 
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Figure  41 
Blue  Exposure  =  Vassec;  Green  Exposure  =  Vioosec ; 
initial  temp  =  -10°C 


Figure  42 
Blue  Exposure  =  6sec;  Green  Exposure 
initial  temp  =  -10°C 


=  2sec; 
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Figure    43 

Field  Power  vs.  Corner  Distance  CX  ; 
r»-.  4^-45°    \  =  20cm    Y,  =  K,X  =  5c 
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Figure   44 
Field    Power    vs.    Corner   Distance   CX; 
y'=-.9Z-30o  >.=    20cm         Y      =•   K,  X  =    5cra 
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Figure    45 
Field    Power    vs.    Probe   Distance   K,X: 
p  =    -.  4^-45°  "X=    20cm         D    =    C*=    10 
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Figure   46    . 
Fi^ld    Power    vs.    Probe    Distance    K.> -  . 
?=   -.9^-30°         "X=    20cm        D   =    CX  =    10 
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coefficient   (p  )   equal   to   -.4Z-45*.   Figure  44  shows  the 
same  for  a  reflection  coefficient  of  -,9Z-30°. 

Figure  45  and  46  show  the  plots  for  the  same  reflection 
coefficients  respectively  as  a  function  of  probe- to- driven 
element  distance.  These  show  the  importance  of  providing 
good  electrical  connection  between  reflector  and  ground 
plane. 
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VII.   EXPERIMENTAL_CONCLUSIONS 

A  particular  set  af  parameters  was  selected  as  being 
optimal.  The  set  and  the  bases  for  its  selection  are 
presented  below. 

A.  FILM  SENSITIZING  PARAMETERS 

1 •   Light  Source  and  Intensity 

The  choice  of  a  4800°K,  225  foot-candle  source  was 
principally  an  economic  decision  based  on  available  devices 
and  a  consideration  of  the  discussion  in  appendix  A. 

2.  Filters 

The  choice  of  Wratten  47B  (Blue)  and  Hratten  99 
(Green)  filters  was  based  on  the  corresponding  responses 
indicated  in  Figures  8  and  9  and  on  the  accompanying 
discussion  in  Section  IV  D.  These  choices  were  considered 
to  result  in  approximately  equal  sensitization  of  blue  and 
green  dyes. 

3.  Sen si tizin^_ Pre- exposure 

Blue  and  Green  pre-exposures  of  2  and  1  second 
respectively  were  chosen  on  the  basis  of  the  discussion  in 
Section  VI.  D  to  place  the  response  in  an  easily 
controllable  region  of  the  surface  of  Figures  37  and  38. 
The  shutter  opening  was  f/5.6  with  a  50mm  lens  focused  to 
60  feet.   The  camera  was  a  Leicaflex. 

B.  FILM  RADIATION  PARAMETERS 
1 .   Initial  Temperature 

An  initial  temperature  of  -20°C  was  chosen  to 
reduce  the  effect  of  ambient  temperatures  while  reducing  the 
amount  of  titre  needed  to  be  spent  preparing  the  film.  As 
indicated  by  Figure  33,  this  temperature  is  considered  to 
have  been  reached  after  1m  and  45s  of  cooling.  Figure  35 
indicates  that  after  30s  of  exposure  to  ambient  temperatures 
the  film  would  still  be  at  less  than  5°C/  and  after  45* 
would  still  be  at  less  than  10°C. 
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2.  Kadiat ion  Frequency 

A  frequency  of  1.5  GHz  was  chosen  to  permit  the 
largest  possible  image  for  the  desired  target  elements  on  a 
single  film. 

3.  Fa diat ion _  Exposure 

It  is  concluded  that  ultimately  the  process  is 
power  limited  with  the  facilities  available  to  this  author. 
At  no  time  was  there  available  a  power  density  greater  than 
approximately  300mw/cm2.  With  that  limitation,  the  only  way 
to  increase  applied  energy  was  to  increase  radiation  time. 
This  was  counter-productive,  since  increased  radiation  time 
allowed  increased  warming,  which  resulted  in  images  being 
fogged  with  responses  to  warmer  temperatures.  A  choice  of 
45  seconds  of  radiation  with  the  maximum  power  available  of 
approximately  16  watts  was  made  to  optimize  the  process  in 
consideration  of  these  factors. 

C.    ALTERNATIVE  SOLUTIONS 

Within  the  overall  aim  of  providing  a  simple  and 
relatively  inexpensive  method  of  imaging  fields,  the  most 
obvious  method  to  improve  the  images  beyond  what  vas 
accomplished  here  would  be  to  increase  the  available  power 
level,  which  would  increase  the  ratio  of  desired  to 
undesired  thermal  fields  in  which  the  process  operates.  By 
thus  raising  the  signal,  the  signal  to  noise  rate  is 
improved. 

The  next  more  obvious  alternative  is  to  reduce  the 
noise  by  reducing  the  ambient  temperature.  This  would 
require  conducting  the  experiments  in  a  cooled  environment 
large  enough  to  contain  the  radiation  structure,  or  at  least 
the  scatterer  and  film. 
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VIII 


FINAL  RESULTS 


A  series  of  images  were  obtained  using  the  parameters 
of  Chapter  VII.  The  target  scatterers  correspond  to  devices 
measured  by  Eurton  [2].   They  are  presented  below. 


//// 1 

>  t  n  ' 

%\ 


Figure  47 
Eegion  of  Cross  Included  in  DTM  Images 
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Figure  48 
DTM  Image  of  Driven  Cross 
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Figure  49 
DTM  Image  of  Scattering  Cross  Y1  =  10cm 
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IX.   CONCLUSION 

A.    SUMMARY 

Dye  Transfer  Modulation  has  been  analyzed  and 
parameterized.  An  optimum  set  of  parameters  subject  to 
certain  constraints  has  been  identified.  Results  have  been 
obtained  on  significant  structures  using  those  parameters. 

A  single  major  problem  remains  to  be  solved  for  the 
process  to  have  widespread  application.  That  problem  is 
that  of  insufficient  margin  between  power  density  of  the 
applied  electromagnetic  field  and  power  density  of  the 
ambient  thermal  field. 

Two  solutions  have  been  suggested,  one  being  to 
increase  the  available  EM  power  density,  the  other  being  to 
decrease  the  ambient  thermal  field  by  conducting  the  process 
in  a  cooled  enclosure. 

A  third  alternative  was  determined  shortly  before  the 
completion  of  this  work.  A  long  wavelength  (2  -  5.6arn) 
infrared  detector  was  obtained  on  a  loan  basis  from  the 
manufacturer.40  The  device  is  sensitive  to  temperatures  in 
the  ambient  range  and  displays  them  on  a  television  screen. 
An  experiment  was  conducted  by  the  author  to  determine  its 
applicability  to  our  purpose.  A  piece  of  ordinary  carbon 
paper  was  attached  to  a  cardboard  backing  with  adhesive.  It 
was  then  placed  against  a  raonopole  scatterer  which  was  then 
irradiated.  The  image  produced  on  the  screen  was 
photographed.   Figure  50  is  the  result. 

The  advantage  lies  in  the  fact  that  it  expresses  a 
change  in  energy  as  a  change  in  hue  of  the  image  rather  than 
simply  a  change  in  density.  Since  the  eye  is  more  responsive 
to   changes   in   color   than   in   density  the  image  presents 


♦©The  device  is  the  AGA  Thermovision  470  System, 
manufactured  by  AGA  Corporation,  550  County  Ave.,  Secausus, 
N.J.,  07904. 
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greater  visibility.   Comparison  between  Figure  50  and  Figure 
49  are  invited. 


figure  50 
AGA  Thermovision  Image  of  Scattering  Cross 

Follow   up   experiments   were  conducted  with  this  device  and 
are  discussed  in  in  Appendix  F. 
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APPENDIX_A 
BASIC  CONCEPTS  RELEVANT  TO  FILM  SENSITOMETRY 

A.    FILM  SENSITIVITY  VERSUS  EXPOSURE 

Direct    parallels   can   be   drawn   between  fundamental 

radiometric  and   photometric   quantities.   In   sensitometry , 

("the    science   of    measurement   of   the   sensitivity   of 

photographic   materials")*1     quantities    which    are   of 

particular  interest  are*2 

Quantity  Symbol  Defining  ..equation   Representative  Units 
jTarae  ~~  "    Ra.2i£i!i^l£ic  ll£l2^etric 

Energy     CX  joule       lumen-second 

(Ira.  sec) 
Energy.    w       u>««*Q/a„      joule/m3     lra.sec/m3 
density  ' °v 

Energy      ^      $.  <*Q/dt       watt        Ira. 

Irradiance  _.       __  watt/m2     lm/m2=lux (lx) 

(illumi-)  t.       t-  ^$/M  =foot-candles 

nance) 

Intensity   I       X'  a$/£w       watt/ste-   candela  (cd) 

ult,^..  j  *„.,  i,  radian     =lm/s teradian 

w*soi.a»^it      (w/str) 

Radiance    ,  .  watt  per     cd/m2 

(Lumi-    i-       L-oX/^w       (strea2) 
nance) 

The  distinction  between  the  two  measuring  concepts  is 
simply  the  preposition  of  the  appropriate  adjective;  e.g. 
radiant  energy  in  radiometry,  luminous  energy  in  photometry. 
This  method  of  distinction  is  evidenced  in  the  quantities  E 
and  L  above.  The  candela  was  defined43  by  the  9th  General 
Conference  on  Weight  and  Measures,  in  1948,  as  1/60  the 
luminance  of  a  black  body  radiator**  heated  to  the 
temperature  of  solidification  of  platinum.  A  lumen   is  the 

♦*Engel,  p. 21. 

♦2"Section  2  -  Radiometry  and  Photometry"  -£.£§£ 
Handbook  of  Photoqraphic  Science  and  Engineering,  1973, 
p  p7~T3o=  1 4BT 

*3Engel,  p.  178. 

♦♦The  subject  of  Densitometry  is  discussed  at  length  in 
"Section  15,  Densitometry , "  J.  Paul  Weiss,  editor,  SPSE 
Handbook  of  Photographic  Science  and  Engineering,  Wiley, 
T9717~pT  8-297~H" 
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time  rate  of  energy  radiation  (radiation  power  output 
present)  in  one  steradian  from  a  source  with  an  intensity  of 
one  candela.  Illuminance  (or,  illumination)  is  commonly 
measured  in  foot-candles  in  the  English  speaking  world, 
which  is  defined  as  one  lumen/ft2.*5  Illuminance  varies  as 
the  square  of  the  distance  from  the  source  to  the  film.*6 
One  more  quantity  is  the  quantity  exposure  (H)  .  Exposure  is 
defined  as*7  an  energy  density: 

H  =  E  t 
where  E  is  defined  above  and  t  is  time.  Exposure  is 
conventionally  given  in  units  of  ergs  per  sguare  centimeter. 
In  practice  it  is  referred  to  as  a  combination  of  lens 
setting  and  exposure  time  in  fractions  of  a  second,  e.g. 
"f/8  at  1/30."  Lens  settings  are  measured  in  "f-stops"  where 
the  f-stop  number  is  defined  as  the  ratio  of  focal  length  of 
the  lens  to  the  diameter  of  the  entrance  pupil.*8  F-stops 
are  conventionally  1.4,  2.3,  4,  5.6,  8,  11,  16,  22  with  one 
half  stop  in  between.  While  most  cameras  are  provided  with 
"click  stops"  which  provide  specific  stop  settings  some 
lenses  may  be  varied  continuously.  One  conventional  f-stop 
decrease  in  lens  opening  represents  an  approximate  doubling 
of  illumination,  hence  of  exposure.  Precisely,  the 
comparison  is  between  the  squares  of  the  f/number:*9 

For  example:  J 

f/16  to  f/11  yields  It=xJ16]2  =  2.121, 

Over  a  short  range  of  f-stops  the  increase  of  an 
additional  f-stop  may  be  compensated  for  by  cutting  the 
exposure  time  in  half  providing  a  constant  exposure  (H)  . 
However,  if  a  large  change  in  either  parameter  is   made   it 

*SEngel,  p.  178-179. 

*6I^id.  p.  79. 

♦7SPSE  Ha nd bq ok,.^ ,  p.  767. 

*8Ibid.  p.  195. 

**Engel,  p.  96. 
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cannot  be  so  simply  compensated  by  changing  the  other 
parameter  by  a  reciprocal  amount.  The  reason  for  this 
failure  of  reciprocity  lies  in  the  reaction  of  the  emulsion 
to  very  large  or  very  short  durations  of  illumination,  with 
different  explanations  for  each  case.  In  either  case  the 
phenomenon  is  referred  to  as  the  Reciprocity  Law  Failure. 5° 

Film  sensitivities  are  conventionally  provided  as  a 
plot  of  experimental  data: 

The  particular  emulsion  type  is  subjected  to   different 

amounts    of   exposure   (H)   and   developed   under   standard 

conditions.  The  density (D)  of  the  developed  negative  is   the 

logarithm   of  the  opacity  (0)  where  the  opacity  is  the  ratio 

of  incident  intensity  to  transmitted  intensity: 

0  =   I  i-Il£i2.~nt_   ,  D  =  log   (0) 
T  transmuted  l  ° 

In   practice,  Density  can  be  quantified  by  visual  comparison 

of  a  projected  negative  image  to  a  card  standard   which   has 

been   printed   with  different  tones  of  grey  corresponding  to 

different  densities.51 

A     plot   of    density   versus   log   H   is   then   the 

conventional  method  of  documenting  the  sensitivity  of  a  film 

to   different   intensities.  Figure  51  gives  a  characteristic 

•'Density  versus  log  H  curve"  (H  and  D   curve)  ,   showing   the 

standard   features.    Various  classification   methods   have 

chosen  cne  or  more  of  the  features  of  the  D-log  H  curve  as  a 

criteria  for  specifing  some  index  to  the  film's  sensitivity. 

In  this  country  the  most  common  is  the   "ASA"   rating   which 

has   been   established   by   the   American  National  Standards 

Institute  (ANSI)  and  generally  is  given  by 

K 

ASA  Speed  =  H 
m 

K  is  a  constant  and  H   is  a  particular  value  of  H  which  cor- 

m 


soibid. ,  pp.  97-93. 

51The  subject  of  Densitometry  is  discussed  at  lenath  in 
"Section  15,  Densitometry,"  J.  Paul  Weiss,  editor,  . SPS E 
Handbook  of  Photographic  Science  and  Engineering,  Wiley, 
797T7"p7  829,  If: 
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relates   to   desired   results  for  different  emulsion  types. 

The  H  and  D  curve  for  a  representative  color   reversal   film 

is  given  in  Figure  52.   For  this  emulsion  type, 

H   =  /H  H 
m   V  1  s 

and 

ASA  Speed  =  H 
m 

rounded  to  the  nearest-3/?  step.52 

B.    FILM  SENSITIVITY  VERSUS  COLOR  OF  LIGHT 

A  different  parameter  of  film  sensitivity  is  the 
frequency  response  of  the  film  emulsion.  Color  sensitcmetry 
provides  curves  of  spectral  density  versus  wavelength  for 
the  individual  layers  of  the  film,  or  alternatively,  an 
impressive  array  of  different  methods  of  quanitizing  the 
response  of  a  particular  film.53  Of  particular  interest  is 
the  concept  of  color  balance.  As  a  rule,  one  desires  a  film 
to  reproduce  the  color  of  an  object  as  it  was  seen.  But  for 
reflecting  objects  the  color  seen  depends  on  the  color  of 
the  illuminating  source.  Therefore,  films  are  constructed  to 
provide  a  "natural"  response  to  objects  illuminated  by 
source  of  a  particular  color  temperature.  Natural  response 
is  considered  to  exist  when  neutrals  are  reproduced  as 
nearly  neutral.5*  The  color  temperature  of  the  source  is 
the  temperature  at  which  a  black  body  radiator  would  emit  a 
spectrum  of  radiation  similar  to  the  source  in  question.55 
If  a  film  is  exposed  to  light  from  a  source  of  a  color 
temperature  other  than  that  for  which  the  film  has  been 
designed  colors  will  not  be  rendered  natural  in  the 
resulting  image.  Thus  the  familar  yellow  tinge  to 
photograghs   taken   under  tungsten  lighting  with  "Daylight" 

S2SPSE  Ha ndb_o_c>k._.ii ,  pp.  810  -  816. 

"Ibid.,  P-  4^8  -  469. 

5*Ibid.,  p.  469. 

55Eastman  Kodak  Company,  Kodak  Publication  No.  U-72, 
Kodak  IFTRAN  Infrared  0£tical  listeria  ls^  1971,  p.  5. 
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Figure   51 
Typical    H   and   D   Curve 


*"*—   (-l'^k#n-   Point  o^  TdnafncV 


Ci/rvC    "Tan^en-t-    T^rcj,,  ^       Point  L. 


L  "*2 —  Berisrt«A    0.£   abolje   Mmimurd  D^ify 


MinimuMD^"^*^ 


Figure    52 
ASA  Speed  Parameters 
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film.  This  imbalance  can  be  corrected  through  the  use  of 
filters  whose  spectral  transmittance  is  such  that  the 
spectral  intensities  of  light  transmitted  through  the  filter 
corresponds  to  the  spectral  sensitivity  of  the  film  in   use. 

C.    FILIEBS 

Several  classes  of  optical  filters  exist  of  which  the 
class  absorption  filters  are  of  interest  here.  Using  the 
guantity  flux  (§)  defined  above  various  gualities  of 
absorption  filters  can  be  considered:56 

Spectral  Reflectance  p-x  is  the  ratio  of  flux  reflected 
by  a  material  (sr)  to  flux  incident  on  it  (s^)  : 

Spectral  Absorption  ft>  is  the  ratio  of  flux  absorbed  by 
the  material  ($$,)  to  the  incident  flux: 

«u  ■  *•/«< 

Spectral  Transiiiittance  y^  is  the  ratio  of  flux 
transmitted  through  the  object  (£c)  to  the  incident  flux: 

"3^.  max  is  the  maximum  transmittance. 

Spectral  density  D>.  is  defined  as 

The  spectral  density  of  two  or  more  filters  used 
simultaneously  is  equal  to  the  sum  of  the  individual 
spectral  densities.57  The  terra  "Spectral..."  implies  a 
dependence  of  the  guantity  on  the  wavelength  of  the  energy 
being  considered.  The  bandwidth  or  half- width  a\,.  is  defined 
as  the  wavelength  interval  between  the  wavelengths  at  which 
T-x  falls  to  iT>  max  and  mean  wavelength  "X^  is  the  wavelength 
about  which  the  pass  band  is  centered.  Of  the  absorption 
filters,  two  types  are  of  interest,  the  band-pass  filters 
and  the  colcr  conversion  filters. 


S6S?3S  Handbook ...  f    p.  259 
57Ibid.,  p.  302-303. 
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1 .   Band-pass  filters 

The  more  common  band-pass  filters  are  constructed 
of  a  gelatin  film  or  colored  glass.  Those  produced  by  the 
Eastman  Kodak  Company  are  classified  by  a  number  series,  the 
"Wratten"  numbers.  These  filters  permit  transmission  of  a 
band  of  wavelengths,  and  may  be  tabulated58  by  mean 
wavelength  X>^  ,  The  Wratten  numbers  increase  approximately 
with  increasing  X^.  Figure  9  above  shows  representative 
half-widths  for  selected  filters.  An  important  aspect  of  a 
particular  filter  is  its  peak  transmittance  (O^max) . 
Different  filtering  materials  have  different  spectral 
transmittances  at  different  frequencies.  Thus  one  obtains  a 
transmittance  curve  for  a  given  filter  similar  to  Figure  53. 
The  differences  in  peak  transmittance  among  the  various 
filters  gives  rise  to  the  requirement  that  different  times 
be  used  to  obtain  equal  negative  density  (D)  when  using 
filters  with  different  peak  transmittances.  For  example,59 
a   Wratten   18A   filter   has   a  >>  max   of   0.62    while    a 

T4 


OA» 


OA 


OZ- 


> 


4.00 


6oo 


too    n.m. 


Figure  53 
Transmittance  of  a  Wratten  47  Filter 


seibid.,  p.  308-309, 
S9ibid. 
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Wratten  99  filter  has  a  X  max  of  0.20  (approximately).  Then 
to  obtain  an  equal  exposure  at  their  respective  should 
require  an  exposure  time  three  times  as  long  with  the 
Wratten  99  filter  as  with  the  Wratten  18A.  The  approximate 
values  of  "D>  max  for  Wratten  47B  and  61  filters  are 
respectively,  .50  and  .40  indicating  exposure  times  of  1.25 
longer  for  the  green  (Wratten  61)  filter. 
2 •   Colcr  Conversion  Filters 

The  need  for  matching  color  balances  between  film 
and  source  tc  render  "natural"  colors  in  the  image  is 
satisfied  through  use  of  color  conversion  filters,  which 
adjust  intensities  of  incoming  light  to  provide  a  wavelength 
intensity  distribution  which  matches  the  sensitivity 
distribution  of  the  film  in  such  a  way  as  to  cause  the 
source-film  system  to  have  an  overall  sensitivity  similar  to 
that  of  the  eye.60  There  exists  a  series  of  filters 
designed  to  provide  this  spectral  shift.  Each  of  the  series 
is  rated  by  its  mired  shift  value.  This  value  (tiSV)  is 
given  by61 

MSV  =  10&  _  io6 

2         1 

where  T   is  the  color   temperature  of  the   source   and  T   is 

1  2 

color  temperature  of  the  light  transmitted  by  the  filter 
(HIRED  is  an  acronym  for  million-reciprocal-degrees)  .  For 
example,  if  a  film  which  has  been  balanced  for  5600°K  is 
used  with  a  4800°K  source  a  mired  shift  of 

10&  _   1 c&  =  -29.76 

is  indicated.  This  corresponds  roughly  to  a  Wratten  82B 
filter,  whose  Mired  Shift  Value  is  given  as  -32. 62  Thus  an 
increase  in  temperature  requires  a  decrease  in  Mired  Value. 


60£ngel,  p.  44. 

61SPSE  Handbook ,    p.  312-317, 

62Ibicl.,  p.  315. 
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POLAROID  CORPORATION 

CAMBRIDGE,  MASSACHUSETTS  02139 


Thank  you  for  calling  us  on  the  Polaroid  Technical  Assistance 
telephone  line. 

Enclosed  is  the  literature  we  discussed.   I  trust  that  you  will 
find  it  informative.   If  you  have  any  questions  in  this  area,  or 
if  you  would  like  help  with  any  other  Polaroid  industrial  product, 
film  or  application,  please  call  us  again  on  our  "Hot  Line." 

We  appreciate  your  interest  in  Polaroid  photographic  products. 

Sincerely, 


In  response  to  your  questions  regarding  the  slope  and  re- 
solution  of  Polaeolor  film,  I  am  pleased  to  pass  on  the  follovin, 
information: 

Slop_e  (of  the  neutral  column  measured  between  .3  and  .7 
densities  using  #92,  93  and  94  Wratten  Statu.;  A 
Filters) 

Red  -  1.32  ±  .15 
Green  -  1.40  ±  .15 
Blue  -     1.24  ±  .15 

Resolution  (contact  print)  -  13  line  pairs/; 


mm 


I  hope  that  this  data  along  with  the  spectral  sensitivity  in- 
formation given  to  you  over  the  phone  will  prove  helpful  in  your 
work .  u 


Fhe  attached  K&D  curve  i s _r o rjiCLagrLtii^l_vp  of  the  Type  10S 
Polaeolor  film  as  measured  by  the  Polaroid  system.   The  initial 
step  on  the  horizontal  exposure  axis  has  an  exposure  of  tv;o 
meter  candle  seconds  from  a  xenon  strobe  light  source  filtered 
to  represent  average  outdoor  daylight.   The  exposure  scale  is 
in  density  units  (logarithmic  scale)  so  the  exposure  at  1.0 
density  would  be  0.2  m.c.s.   The  vertical  scale  is  the  sample 
density  as  measured  on  a  reflection  densitometer. 

The  spectral  sensitivity  data  attached  is  from  a  studv  done  in 

our  laboratories  several  years  a<~o  '  us  In?;  a  non-equal  energy 

spectrophotometer  on  separate  film.. samples  for  each  wavelength 

measurement .   The  data  is  plotted  as  a  function  of  wavelength 

against  the  I03  of  the  number  of  incident  photons  cor   square 

centimeter  to  reach:  a  1.0  density  on  the  positive.  -.Although  v:e 

do  not  have  an  estimate  of  the  error  involved  in  this  technique, 

the  data  is  a  reasonable  representation.   Subsequent  to  this 

studv,  a  modification  of  the  nerative  has  been  introduced  into 

the  product  to  reduce  the  green  sensitivity  in  the  UV  and  blue 

region  of  the  soectrum.        •'  ,       •   ■ 

*     •    •  ... 

both  pieces  of  attached. data  are  from  single  film  lots  represen- 
tative of  the  Type  103  Polaeolor  product.   As  such.,  they  nho-:':^ 
not  be  r  -i  scons  trues'  as.  an  avera  ~y  e  r  e  s  u  1 1  r.  c  r  as  ,a  product 
Kneel  ??.  ca  t  i  o n  .       .         "     ~ 
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APPENDIX. C 

COMPUTER_PROGRAM 

A.  DISCUSSION 

Program  PWRCK  computes  the  values  of  power  versus 
corner  to  monopole  (CX)  and  monopole  to  probe  (K,X) 
distances  as  described  in  section  V.  C.  2.  The  data  thus 
computed  is  displayed  as  either  a  two-  or  three-  dimensional 
plot.  The  type  of  plot  generated  depends  on  the  parameters 
IC  and  IK.  If  either  of  these  is  set  equal  to  the  integer  1, 
a  2-D  plot  shall  be  generated.  Otherwise  a  3-D  plot  is 
generated. 

Since  the  purpose  is  to  provide  flexibility  in 
computation  there  is  necessarily  a  finite  complexity  to  the 
program.  This  primarily  comprises  determining  the  various 
parameters  of  interest  and  inserting  them  in  the  program. 
Insertion,  rather  than  reading,  of  the  values  has  been 
selected  for  two  reasons.  First,  computer  time  is  reduced, 
and  second,  once  set,  the  majority  of  the  parameters  are  not 
likely  to  be  changed. 

Only  one  other  user  interaction  is  required  besides 
setting  parameter  values.  That  is  the  task  (unavoidable  in 
FORTRAN)  of  explicitly  defining  the  size  of  the  arrays. 
Three  sets  of  arrays  must  be  DIMENSIONed.  The  first  deals 
with  the  power  computations.  The  second  deals  with  plot 
scaling  for  three  dimensional  plots.  The  third  set  affects 
the  output  routines  but  is  fixed  and  does  not  require  user 
interaction.  Procedures  are  described  below. 

B.  PARAMETERS  TO  BE  PROVIDED 

The  following  values  are  required  for  all  types  of 
output: 

CAPH  =  CAP  H  =  Monopole  length  PI  in  centimeters 

CLMDA   =   Lambda   -  Radiation  frequency's  wavelength  (>) 
in  centimeters 
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1C  -  number  of  data  points  computed  for  each  P  versus  C 
curve 

IK  =  number  of  data  points  computed  for  each  P  versus  K 
curve 

12  =  number  of  integration  steps  per  datum  point 

RC  =  magnitude  (always  positive)  of  the  complex 
reflection  coefficient  P. 

ARC  =  phase  angle  of    in  degrees 

Y1  =  initial  monopole  to  probe  distance  K  x  in 
centimeters  i 

YMAX  =  largest  value  of  K,X  in  centimeters  by  which 
data  is  to  be  computed 

DY  =  Interval  in  K,\  in  centimeters  by  which  data 
points  are  to  be  separated 

EZO  =  an  arbitrary  scale  value 

D  =  Initial  corner  to  monopole  distance  CX  in 
centimeters 

DMAX  =  largest  value  of  CX  in  centimeters  for  which 
data  is  to  be  computed 

DD  =  Interval  in  CX  in  centimeters  by  which  data 
points  are  to  be  separated 

Z  =  Probe  length  (Sx)  in  centimeters 

DZ  =  Integration  interval  in  Z  with  which  Simpson's 
rule  integration  is  performed 

TTL  =  a  96  character  alpha-numeric  title  to  be  read  in 
on  two  successive  cards,  with  48  characters  (blanks 
included)  per  card 

C.    DETERMINATION  OF  PARAMETERS 

1  •   Defined. parameters 

Certain  parameters  are  defined  in  metric  units 

CAPH 

CLMDA 

RC 

ARC 

Y1 

EZO 

D 

YMAX 

DMAX 

Z 

Note  that  YMAX  and  DMAX  are  not  explicitly  required   by   the 

program. 

2-   Derived  Parameters 

In  all  of  the  derivations  the  issue  is  to  determine 

the  interval  between  computations  in  each  dimension  (DY,  DD, 

DZ)   which   gives   the   smoothest  results.  The  computer  must 
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also  be  told  the  number  of  such  steps  (IK,  IC,  L2)  which  are 

being   taken.   There  is  no  way  to  avoid  mental  gymnastics  in 

deriving  these  values  since  one  must  be   determined   on   the 

basis   of  the  other  in  an  iterative  manner.  It  is  considered 

that  this  is  most  efficiently  done   off-line,   as   described 

with  the  below  formulas.  It  is  worthwhile  to  note  that  there 

is  always  one  more  in  the  number  of  data  points  than  in   the 

number   of   intervals  in  a  range,  since  one  must  account  for 

the  very  first  (zero-interval)  computation.   On  these   bases 

the    following    formulas   generate   the   required   derived 

parameters. 

IC  =  DMAX  -  1  +1  DD  =  DM AX  -  DP 

DD  IC  ^~~T 

IK  =  YMAX  -  Y1  +1  DY  =  DM AX  -  DP 

D¥  IK  -  1 

L2  =  Z_  +  1  DZ  =    Z 

DZ  L2  -  1 

D.    SAMELE  PROBLEMS 

The   three   sample   problems  demonstrate  the  use  of  the 
program. 

1  •   Two _  Dimen s ig na  1^ P 1  ot (Power  ys_Y(-p_f prefixed _D) 

Consider  a  problem  in  which  a  two  dimensional  plot 
is  desired  cf  Power  versus  Y,  for  a  fixed  value  of  CX  with 
a  given  reflection  coefficient  -.4^-45°.  Refer  to  Figures 
17  and  18  as  desired.  A  wavelength  of  20cm  (corresponding 
to  a  freguency  of  1.5  &Wz)  is  selected  with  a  fixed 
Corner-to-  Monopole  distance  of  5cm,  probe  length  of  2cm  and 
raonopole  length  (V4)  °f  5cm.  Initial  value  of  Y,  is  2cm 
and  data  is  desired  out  to  one  wavelength  away  from  the 
monopole. 

a.  Defined  parameters  from  this  problem  statement 
are 

CAPH  =  5.0 
CLMDA  =  20. 
RC  =  0.4 
ARC  =  135. 
Y1  =  2. 
EZO  =  30. 
D  =  5. 
YMAX  =  20. 
DMAX  =  5. 
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Z  =  2. 
b.   Derived  Parameters 

Only  a  single  curve  is  desired,  and  there  is  a 
single  value  of  D,  thus  IC  =  1.  In  order  to  get  a  smooth 
plot  it  is  likely  that  an  interval  of  about  20  points  per 
inch  should  give  a  smooth  curve.  But  to  allow  easy 
graduation  of  the  axes  one  ought  to  have  some  metric 
relationship  between  intervals  and  the  range  YMAX  -  Y.  Since 
in  this  case  YMAX  -  Y  is  18cm.,  a  choice  of  90  intervals  is 
suggested,  giving 

90  =  YMAX  -  Y  so  that  IK  =  91 
DY 

Similarly,  a  choice  of  approximately  100   integration   steps 
results  in  L2  =  10  1.  From  these  we  obtain 


IC  =  1 

IK  =  91 

L2  =  101 

DD  =  1.0 

DY  =  0.2 

DZ  ^  0.02 

With  these  parameters  inserted  and   the   title   cards   (TTL) 

appended   in   the   conventional   manner   as   data  cards,  the 

listing  on  page  105  results  in  the  5"  x  3"   plot   of   Figure 

54. 
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Figure   54 
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2.   Two  Dimensipna l_Plot  XPQwer_yersus_ D  f or _f ixed  Y  ) 
a.   Defined   Parameters 

In  this  case,  using  the  same   scenario   as   (1) 
above,  but  plotting  versus  the  other  parameter  requires  only 
the  change  in  Y1   and   D.   For   example   if   a   fixed   probe 


distance   Y1 


o. 


18cm   were  selected  and  the  initial  corner 


distance  D  were  selected  to  be  2cm,  the   two   cards   changed 
would  be 

Y1  =  18.0 

D  =  2.0 

b.   Derived   parameters 

Since  the  intervals  involved  have  not  changed, 
that  is,  probelength  is  the  same,  and  the  data  line  is  still 
18cm.  The  computations  are  conceptually  identical  with  the 
first  problem.  Axes  have  changed,  so  the  intervals  and 
counters  change 

'  IC  =  91 

IK  =  2 
L2  =  101 
DD  =  0.2 
DY  =  1.0 
DZ  =  0.2 

The  resulting  plot  is  shown  in  Figure  55. 
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Figure    55 


100 


3  •   Ik£§e_J)imensional_Plot._XPow  er_  vers  us_Yj a n d_D^ 

a.  Combining  the  reasoning  of  the  first  two 
problems  one  might  choose 

Y1  =  2.0 
D  =  2.0 

b.  Derived  Parameters 

Now,  unfortunately  the  number  of  computations 
has  been  squared  unless  some  adjustment  is  made.  The  three 
dimensional  plot  is  desireable  primarily  as  an  aid  to 
visualization  rather  than  for  scientific  accountability,  so 
it  is  reasoned  that  a  less  precise  plot  is  acceptable.  Here 
two  intervals  per  centimeter  of  data  range  is  considered 
adequate.  Similarly,  the  number  of  integration  steps  are 
reduced,  so  that 

IC  =  37 

IK  =  37 

L2  =  51 

DD  =  0.5 

DY  =  0.5 

DZ  =  0.04 

c.  For  the  three  dimensional  plot,  some 
additional  user  interaction  is  required  by  the  plotting 
subroutine.  These  are  described  in  detail  in  the  Technical 
Memorandum  for  PLD3D1  available  from  the  Computer  Center 
Consultant  Staff.  This  program  has  reduced  the  parameters 
thereof  tc  three  items. 

NKXY  -  a  "magic  number"  which  is  determined  by  PLT3D1 

ALPHA  -  rotation   angle   of   the   3-D   surface   about   a 
horizontal  axis 

BETA  -  rotation  angle  of  the  surface   about   a   vertical 
axis 

These  three  values  are  to  be  read  in  when   provided   in   the 
above  order  in  an  (110,  2F10.3)  Fortran  FORMAT. 

(1)  Initially  NKXY  is  set  to  -1  and  entered 
in  column  9.  ALPHA  and  BETA  must  be  set  to  their  desired 
value.  The  program  provides  two  output  plots,  so  two  data 
cards  must  be  provided.  It  is  suggested  that  the  angles  be 
selected  initially  to  be  15°  and  15°  respectively  for  ALPHA, 
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15°   and   75°   respectively   for   BETA.  Accordingly  the  data 
cards  should  appear  as 

- i      15.0      i  5 . 0 
l=    J.---  i  i\Aiii_-irrii_  MiM •!■•  ?  KL— «  4 ■•'■. I J5 

-i      15.0      P5.0 
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8  0  0 ,"  P  0  3  0  0  0  0  0  D  :  8  0  0 ,' ,-  G ,?  P  C  3  0  D  B  C  0  0  ■  •  0  0  0  C  0  0  B  0  0  0  0  0 

'    :    J    4    s    '    '    !    J   '«  II   H  II  M  '5  13  IT  13  13  2!  Zl  V.  Z3  24  !S  ZE  :1  II  23  33  Jl  32  J]  )(  »  «  11  !;  )]  II  ||  <•  <] 

M  1!  I  1  1  11  1  1  1  1  1  1  1  II  1  l  1  1 ,;  1  n  1  I  11  1  1  1  1  |  J  l  n  n  1  1 

2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2 ,"  2 ,,  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  I  2  2  2  2 

3  3  3  3  2  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3  3 

4  4  4.  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  ^  4444444444444 

5  5  5  5  5.'  5 ,  5  5  5  5  5  5  6  5  5  5  5  5  5  5  5  5  ~s  55^55^55555555555 
,  G  6  6  6  6  C  6  6  6  S  5j  ,,3  6  6  G  E  6  E  C  C  6  6  6  G  6  6  B,  E  G  G  6  0  G  E  G  6  G  6  E 
7  7,  77777777,  777777777  7777,7777777777777777 
S  8  8  8  8  8  8  8  3  8  8  8  B  8  8  8  8  8  S  8  B  3  8  6  8  8  8  8  „,  8  8 ,  8  8  8  8  S  8  8  8  8  8  8 

9 "  s  s '  s  9  s ;;  5  9  a  9  c  3 "  9  s  s  s  e  s  s  <r  s  3  9  s  ■<  5  s  3  s  s  9  9  9  3  s p  9  s 

'    '    '    <    5    S    1    Z    !   IJ  1112  13  U  15  It  IZ  II  13  25  21  Z2  23  2<  25  25  ."!  28  23  38  21  32  33  31  35  IS  3)  39  23  <0  41  12  <1 
NECC/HP-23431 


(2)   When   ran,   the   program   will   define  the 
required  value  of  NKXY  as  follows 

ALPHA=  15.0      BETA=  15.0 

PLT3D1  WILL   PRODUCE   AN   ISOMETRIC   PROJECTION   OF   THE 
INPUT  SURFACE 

THE  VALUE  OF  NKXY  FOR  THIS  PLOT  SHOULD  BE      0 

SUBROUTINE  PLT3D1  HAS  PLOTTED  A  GRAPH  TITLED: 

***POWER  VS  CORNER  TO  MONOPOLE  DISTANCE  (C)  AND  **** 
***MONOPOLE  TO  PROBE  DISTANCE  (K)-(IN  WAVELENGTHS)  **** 
ALPHA=  15.0      BETA=  75.0 

PLT3D1   WILL   PRODUCE   AN   ISOMETRIC   PROJECTION   OF  THE 
INPUT  SURFACE. 

THE  VALUE  OF  NKXY  FOR  THIS  PLOT  SHOULD  BE      0 

SUBROUTINE  PLT3D1  HAS  PLOTTED  A  GRAPH  TITLED: 

***POKER  VS  CORNER  TO  MONOPOLE  DISTANCE  (C)  AND  **** 
***MONOEOLE  TO  PROBE  DISTANCE  (K)-(IN  WAVELENGTHS)  **** 
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At  this  point  the  user  must  set  NKXY  to  the  proper  values, 
sLUil  check  the  DIMENSION  statement  in  the  procjrara  to  insure 
that  the  arrays  KX  and  KY  are  at  least  as  large  as  NKXY. 

When  these  adjustments  are  made  the  program  is 
re-run,  the  resulting  output  plots  shall  result  as  shown  in 
Figures  56  and  57. 


Figure  57 
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d.   Scaling 

The  program  is  set  up  to  provide  a  real  valued 

output  in  mw  for  an  input  current  Im  of  Im  =  EZO.   It  may  be 

desired  to  normalize  the  plots  to  a   fixed   value   to  allow 

comparison  of  curve  shapes.  This  is  accomplished  by  the 
block  of  source  cards  from  statements  2000  through  2100  as 
fellows: 

2000  CONTINUE 
c . _ C 

c  c 

c  c 

c  c 

c _ _____ _ -_ . c 

DO  2100,  1=1, IC 
P  (N,I)  =  P  (N,I)/PRHX 
2100  CONTINUE 

Other   scaling  techniques  may  be  substituted  for  this  set  of 

statements. 
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APPENDIX_D 
MATHEMATICAL  PROCEDURES 

A.    STATISTICAL   METHODS 

Data  presented  frora  the  various  experiments  are 
identified  as  "average  values"  derived  from  the  stated 
number  of  runs  on  the  experiments.  The  following  algorithms 
were  used  to  obtain  the  various  statistical  values  referred 
to  in  the  text: 

1-   Average  Value  (x)  of  a  set  {x  ,    x  ,...,  x^ }  : 

12 

N 
x  =  £  x: 
i=1 

2.  St andardard,  .Deviation  (Sx)  is  used  to  characterize 
the  accuracy  of  a  particular  average: 


3.  Average^ Standard  Deviation  (S"x)  is  used  to 
characterize  the  overall  accuracy  of  a  curve,  each  point  of 
which  is  an  average  value  with  its  attendant  standard 
deviation.   It  is  the  average  of  those  standard  deviations: 

M 
SX  =  ±    2   SX: 

B.    CURVE  FITTING 

1 .  Temperature_ Curves 

Temperature  Curves  were  fitted  with  a  French  curve 
passing  through  the  average  data  points  for  the  given  set  of 
conditions. 

2.  Density__Curves 

Density  curves  in  Figures  37  and  38  were  adjusted 
within  one  half  a  standard  deviation  on  either  side  of  the 
averaged   datum   point.   These  restrained  adjustments  in  one 
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dimension  were  necessary  to  obtain  consistency  in  the  other 
planar  dimension.  Figure  58  shows  the  technique  used. 
Averaged  data  points  were  plotted  with  a  line  showing  a 
width  equal  to  one  standard  deviation  centered  on  the  datum 
point.  Straight  lines  connecting  these  limits  created  a 
"channel"  (the  hatched  region)  within  which  adjustment  was 
allowed.  Thus  figures  37  through  40  show  values  obtained 
within  one  half  a  standard  deviation  of  the  mean 
experimental  values  with  quasi-linear  interpolation  between 
data  points. 


»en*; 


Fi.tte.cL   Curve 


"Datum 


t>a«i 


Cboat 


Vdtutn 


•°^V\  C^vpofcorc  ) 


Figure  58 
Region  of  Likely  Values  about  a  Datum  Point 
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APPENDIX_E 
EXPERIMENTAL  DATA 
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•  Average  of  3  data  points 
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Figure  37 
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Figure  39 

Initial 

Temp 

Radiation 
Time 

-30.0 

30.0 

-30.0 

32.5 

-30.0 

35.0 

-30.0 

37.5 

-30.0 

40.0 

-30.0 

42.5 

-30.0 

45.0 

-30.0 

47.5 

-30.0 

50.0 

-30.0 

52.5 

-30.0 

55.0 

-30.0 

57.5 

-30.0 

60.0 

-27.5 

30.0 

-27.5 

32.5 

-27.5 

35.0 

-27.5 

37.5 

-27.5 

40.0 

-27.5 

42.5 

-27.5 

45.0 

-27.5 

47.5 

-27.5 

50.0 

-27.5 

52.5 

-27.5 

55.0 

-27.5 

57.5 

-27.5 

60.0 

-25.0 

30.0 

-25.0 

32.5 

-25.0 

35.0 

-25.0 

37.5 

-25.0 

40.0 

Average 
Density 

0. 1565 

0. 1670 

0. 1790 

0. 1905 

0.2000 

0. 2105 

0.2185 

0.2245 

0.2290 

0.2320 

0.2350 

0.2380 

0.2405 

0. 1560 

0. 1670 

0. 1785 

0. 1900 

0.2010 

0.2110 

0. 2200 

0.2  27  5 

0.2350 

0.2425 

0.2465 

0.2520 

0.2550 

0. 1540 

0. 1655 

0. 1780 

0. 1900 

0.2010 


126 


Initial 
Temp 

Radiation 
Time 

-25.0 

42.5 

-25.0 

45.0 

-25.0 

47.5 

t25.0 

50.0 

-25.0 

52.5 

-25.0 

55.0 

-25.0 

57.5 

-25.0 

60.0 

-22.5 

30.0 

-22.5 

32.5 

-22.5 

35.0 

-22.5 

37.5 

-22.5 

40.0 

-22.5 

42.5 

-22.5 

45.0 

-22.5 

47.5 

-22.5 

50.0 

-22.5 

52.5 

-22.5 

55.0 

-22.5 

57.5 

-22.5 

60.0 

-20.0 

30.0 

-20.0 

32.5 

-20.0 

35.0 

-20.0 

37.5 

-20.0 

40.0 

-20.0 

42.5 

-20.0 

45.0 

-20.0 

47.5 

-20.0 

50.0 

-20.0 

52.5 

-20.0 

55.0 

-20.0 

57.5 

-20.0 

60.0 

Average 
Density 

0.2105 

0.2220 

0.2300 

0.2390 

0.2500 

0.2550 

0.2620 

0.2670 

0. 1520 

0. 1650 

0. 1780 

0. 1900 

0.2020 

0.2130 

0. 2245 

0.2335 

0. 2455 

0.2550 

0.2625 

0.2720 

0.2770 

0. 1535 

0. 1660 

0. 1800 

0. 1915 

0.2035 

0.2160 

0.2283 

0.2385 

0.2500 

0.2605 

0.2710 

0.2810 

0.2885 


127 


Initial 

Temp 

Radiation 
Time 

-17.5 

30.0 

-17.5 

32.5 

-17.5 

35.0 

-17.5 

37.5 

-17.5 

40. 0 

-17.5 

42.5 

-17.5 

45.0 

-17.5 

47.5 

-17.5 

50.0 

-17.5 

52.5 

^17.5 

55.0 

-17.5 

57.5 

-17.5 

60.0 

-15.0 

30.0 

-15.0 

32.5 

-15.0 

35.0 

-15.0 

37.5 

-15.0 

40.0 

-15.0 

42.5 

-15.0 

45.0 

-15.0 

47.5 

-15.0 

50.0 

-15.0 

52.5 

-15.0 

55.0 

-15.0 

57.5 

-15.0 

60.0 

-12.5 

30.0 

-12.5 

32.5 

-12.5 

35.0 

-12.5 

37.5 

-12.5 

40.0 

-12.5 

42.5 

-12.5 

45.0 

Average 
Density 

0.  1640 

0.  1750 

0. 1880 

0.2000 

0.2120 

0.2240 

0.2350 

0.2450 

0.2550 

0. 2650 

0.2750 

0.2840 

0.2910 

0. 1760 

0. 1870 

0. 1985 

0.2100 

0. 2220 

0.2330 

0.2440 

0.2535 

0.2630 

0.2720 

0.2805 

0.2890 

0.2960 

0. 1900 

0.2000 

0.2110 

0.2220 

0.2340 

0.2440 

0.2545 


128 


Initial 

Temp 

Radiation 
Time 

-12.5 

47.5 

-12.5 

50.0 

-12.5 

52.5 

-12.5 

55.0 

-12.5 

57.5 

-12.5 

60.0 

-10.0 

30.0 

-10.0 

32.5 

-10.0 

35.0 

-10.0 

37.5 

-10.0 

40. 0 

-10.0 

42.5  • 

-10.0 

45.0 

-10.0 

47.5 

-10.0 

50.0 

-10.0 

52.5 

-10.0 

55.0 

-10.0 

57.5 

-10.0 

60.0 

Averaqe 
Density 

0.2630 

0.2720 

0.2790 

0.2865 

0.2940 

0.2990 

0.2020 

0.2150 

0. 2260 

0.2365 

0.2470 

0.2  57  0 

0.2660 

0.2730 

0.2800 

0.2860 

0. 2910 

0.2950 

0.2990 


129 


-^Figure 

cx^: — — 

43 

—  —  —  ~  T 

6 

51.45 

7 

66.48 

8 

65.86 

9 

63.11 

10 

59.44 

11 

57.56 

12 

54.95 

13 

51.00 

14 

45.55 

15 

41.22 

16 

35.82 

17 

28.95 

18 

23.00 

19 

15.76 

20 

10.72 

5x  =  1.436 

45_ 

10.387 
9.582 

8.739 
8.250 

7.^38 
6.706 
6.116 
5.206 

4.693 
4.586 
4.380 
3.904 
3.643 
3.400 

3.323 


130 


APPENDIX.! 
A GA_THEB MO VISION „ EXPERIMENTS 

A.    DISCUSSION 

The  author  attended  a  general  demonstration  of  the  AGA 
Thermovision  Device  Model  750  on  22  July  1975  at  NPS.  The 
author  recognized  the  possible  application  of  the  device  to 
the  objectives  of  this  work,  and  asked  the  AGA 
representative  to  take  part  in  some  simple  experiments  that 
afternoon  to  see  if  the  application  of  the  device  was 
feasible.  The  representative,  i'ir.  Jack  Patterson,  agreed  to 
do  so.  The  experiment  consisted  simply  of  a  Case  1A  cross 
(described  below)  10cm  from  the  driven  element  irradiated 
with  the  corner  reflector  system  with  the  reflector  1Scm 
from  the  driven  element.  A  Polaroid  Print  backing  cardboard 
on  which  was  attached  a  piece  of  Olde  Towne  carbon  paper  was 
placed  next  to  the  cross.  The  power  density  arriving  at  the 
cross  was  estimated  to  be  approximately  15  raw/cm2  and  the 
ambient  temperature  was  24°C.  The  Thermovision  was  focused 
on  the  cross,  resulting  in  the  image  of  which  Figure  53  is  a 
conventional  Polaroid  photograph. 

Encouraged  by  these  results  and  after  consultations 
with  Professor  Burton,  arrangements  were  made  for  a  more 
extensive  and  controlled  feasibility  test.  Those 
experiments,  conducted  on  17  November  1975,  were  to 
determine  the  feasibility  of  using  the  AGA  Thermovision 
System  to  detect  thermal  fields  induced  by  electromagnetic 
fields  on  an  assortment  of  structures.  The  form  of  these 
experiments  was  selected  to  investigate  two  principal 
elements  considered  likely  to  affect  that  feasibility. 
These  elements  are: 

1 .   Conductivity 

The  surfaces  investigated  must  be  sufficiently 
conductive  to  allow  current  and  charge  distributions  to  be 
established.     In   particular,   current   distributions   are 
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necessary  in  order  to  generate  thermal  fields  which  may  be 
then  sensed  by  the  Thermovision  System.  However/  the 
surfaces  must  not  be  so  conductive  that  either  (a)  no 
significant  heat  is  generated,  or  (b)  any  generated  heat  is 
immediately  dissipated  by  the  thermal  conductivity  of  the 
material.  (Thermal  and  electrical  conductivity  go  virtually 
hand  in  hand  in  electrically  conductive  materials.) 
2 .   Emi  ss  i  vi  t  _y 

The  materials  used  must  have  their  surfaces  treated 
in  such  a  manner  as  to  assure  a  sufficient  level  of 
emissivity. 

a.  Spectral  emissivity  (e%)  is  a  measure  of  the 
difference  between  the  spectral  emittance  of  the  object  in 
guestion  and  that  of  an  ideal  black  body: 

where  W^   is   the   spectral   emittance  of   the   object   in 
question,  Wxs  is  that  of  the  black  body. 

b.  The  emissivity  (g)  or  total  emissivity  of  an 
object  is  the  ratio  of  its  total  emittance  to  that  of  a 
black  body  expected  from  the  Stef an-Boltzraann  formula 

=  0>       __0     =     0 

where  W0  is  the   total   power  density   at   all   frequancies 

emitted   by   the   object,  0"is  the  Stef an-Boltzmann  constant 

(5.67  •  10~8  watt   )   and  T  is  the  temperature  of  the  object 
in  degrees  Kelvin. 

"The  values  for  £  obtained  by  using  Thermovision 
are,  in  effect,  the  average  of  £*  occuring  over  the 
middle  infrared  wavelength  interval  utilized..."63 

The   emissivity   of   a    material   is    highly 
dependent   on   the   condition   of   its  surface.   In  general. 


63Thermovision  Manual,  Section  7,  p.  7. 
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highly  reflective  surfaces  are  very  poorly  emissive,  so  that 
by  increasing  its  absorptance,  we  can  increase  its 
emissivity. 

B.    BASIC  FORM  OF  EXPERIMENTS 

The  experiments  comprised  three  basic  forms. 

1 •   Dielectric  Substrate  .with  Carbon  Paint 

Various  structures  were  made  of  dielectric 
substrate  and  painted  with  a  conductive  carbon  paint.6* 
Carbon  has  been  seen  in  the  first  experiment  to  have  heat 
generating  and  keeping  qualities  which  were  adequate  to 
provide  a  Thermovision  image.  The  question  was  whether  its 
electrical  conductivity  is  sufficiently  high  to  permit 
induction  of  currents  of  sufficient  magnitude  as  to  match 
the  behavior  of  a  more  perfectly  conducting  structure. 

2 .   Me tal_ Structure^ with and  without  Carbon __P a  in t 

These  structures,  matching  those  of  B.1  above, 
were  considered  a  likely  substitute  for  the  problem  of 
inadequate  induced  currents.  The  question  here  is  whether 
they  possess  enough  emissivity  to  present  an  image  to  the 
Thermovision  System. 

3-   Structures  with^Carbon.  Paper 

These  structures  were  similar  to  that  of  the 
earlier  experiment.  The  principle  involved  is  to  use  the 
carbon  paper  to  tranduce  fields  in  the  vicinity  of  the 
structure  from  which  charges  and  currents  on  the  structure 
can  be  deduced.  The  phenomenon  is  considered  to  consist  of 
currents  being  induced  in  the  carbon  paper.  These  currents 
are  parallel  to  the  electric  lines  of  force  (LOF)  set  up  by 
the  charges  on  the  structure.  As  the  lines  of  force 
converge  on  large  charge  concentrations,  their  flux  density 
(lines  per  square  centimeter)  increases,  and  with  it  the 
heat  generated  by  the  ohmic  resistance  of  the  carbon  paper. 
In   high    charge    density    regions    this   LOF    density 


64Eccocoat   256   paint,   manufactured   by   Emerson   and 
Cuming,  Inc.,  Canton,  Massachusettes,  02021. 
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and  its  attendant  thermal  energy  density  (joule/cm2)  is 
sufficient  tc  exceed  the  sensivity  of  the  AGA  detector  and 
provide  an  image. 

C.    EXPERIMENTAL  PARAMETERS 
1.   Types  of  Structures 

Three   basic  structures  were  used. 

a.  1.5  Wavelength  (X)  Monopole  over  a  Ground  Plane 

b.  Aircraft  Model 

c.  Wire  crosses 

2«   Expected _ Charge_and  Current  distributions 

Various  charge  and  current  distribution  are 
possible  on  the  selected  structures.  Broadly  these  may  be 
classed  into  various  cases. 6S 


65Burton/  R. W. ,  "Measured  Currents  and  Charges  on  Thin 
Crossed  Antennas  in  a  Planewave  Field",  (with  R.W.  King) 
IEEE  transactions  on  Antennas  and  Propagation, Vol.  AP-23, 
TTo.""5,  pp.o37-cw,  September  TT7d  . 
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a.   Case  1A 

This  case  uses  a  3/<i"X  vertical  height  (h)  over  a 
ground  plane,  with  cross  arras  (when  present)  placed  V2. 
above  the  ground  plane.  These  arms  are  each  */*♦  to  be 
compatible  with  the  vertical  member's  charge  and  current 
distribution.  The  effect  of  this  configuration  is  to  place 
a  large  current  maximum  at  the  junction  point,  and  large 
charge  deposits  at  the  tips  of  the  structure. 


Vertical 
Current 


Ground 
Plane 


rji ;>  Jr >j,  ?fy 


Nk 


Cross-arm  Current 


V«A 


Vertical 
charge 


Cross-Arm  Charge 

Figure  59  -  Case  1A 
Cross  Configuration  with  Current  and  Charge  Distribution 
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b.   Case  1B 

This  case  also  uses  a  %>.  h  but  places  the 
cross  arm  at  a  charge  maximum.  Now  the  cross  arms  must  each 
be  Vfc  to  have  their  natural  charge  and  current 
distributions  be  compatible  with  that  of  the  vertical 
member. 


9M\ 


■  l  )  i  ijtt  >r  /t  'I 
f  >  I 


Vertical 
Current 


Vertical 
Charge 


Cross  Arm  Charge 

Figure  60  -  Case  1B 
Cross  Configuration  with  Current  and  Charge  Distribution 
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c.   Case   2 

This   case  presents  a  symmetric  structure,  with 
the  vertical  member  one  wavelength  in  length,  and   each  arm 
Va  in   length.    This   results   in  a  current  minimum  at  the 
junction . 


tu 1 1  ni  >  *  > > t 


Cross  Arm  Charge 


Vertical 
Current 


Vertical 
Charge 


Figure  61  -  Case  2 
Cross  Confiauration  with  Current  and  Charge  Distributions 


137 


d.   Case  3 

This  case  presents   a   1.5X  structure   without 
cross  arms. 


i>3/>flttllfttt'l 


'      << 


Vertical 
Current 


Vertical 
Charge 


Figure  62  -  Case  3 
Cross  Configuration  with  Current  and  Charge  Distributions 
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3 .   Dimensions  and_ Frequencies  to  be  used 

The  frequencies  to  be  used  depend  on  the  dimensions 

of  the  structures  involved. 

a.   The   available   materials   for  construction  and 

the  considerations  above  of  the  metal   monopole   required   a 

Case   3   monopole   over   a   ground  plane  with  dimensions  as 

follows: 


a  =  2  inches  =  5.03  cm 


Ka  =  2  a  =  1. 


l*W.ftBcm*I.S\ 


X=  31.92  cm 


f  =  939.39MHz 


1  =  1. 5X  =  47.88  cm 


Figure  63 
Dimensions  of  the  1.5X  Wavelength  Monopole 

From  these  dimensions  an  operating  frequency  of 
939.89  MMz  is  derived. 

b.   Plastic  Substrate  Aircraft  Model 

A  plastic  model  fighter  plane  was  selected  with 
the  dimensions  shown  in  the  figure  below.  From  these 
figures  the  range  of  dimensions  were  obtained  for 


Length  (L) 
Width  (W) 


26.6  4    L  4:   27. 6  cm 
20. 1  *  H  4   25.7  cm 
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Figure  64 
Dimensions  of  the  Plastic  Aircraft  Model 

(1)  A  Case  1A  structure  is  almost  achieved  with 
a  frequency  selection  of  815.2  MHz  with  the  model  nose  down 
over  the  ground  plane.  In  such  a  case  L  =  27.6  cm  =  %"\,  and 
the  arms  are  each  .02   longer  than  XA*  : 


=  », 


ft>  =  2  7.6  cm 


W  =  h\    =  18-4  cm;  20.1  =  .546* 

(2)  A  Case  2  Configuration  is  obtained  at 
1.258  &Mz  with  the  model  tail  down  over  the  ground  plane. 
In  this  case  the  length  is  one  wavelength  but  the  arms  are 
now  each  .02   too  short: 

L  =  *  =  26.6  cm 


W  =  "X  =  26.6  cm; 


25.7  cm  =  .966X 
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c.   Metal  Aircraft  Model 

A  copper  aircraft  model  was  obtained  as  shown 
in  the  figure  below.  This  structure  is  more  symmetric  than 
the  previous  one  so  that  the  only  case  considered  favoraDle 
is  a  Case  2  configuration.   The  ranges  of  dimensions  are 

506  I,   4    54  cm 
50  4:     W  i  57  cm 

A   frequency   of  570  MHz  provides  X=  52.6  cm,  which  matches 
these  dimensions. 


1fek 


Figure  65 
Dimensions  of  the  Copper  Aircraft  Model 

d.   Wire   Crosses 

Wire  crosses  were  constructed  to  provide  Case 
1A  and  1B  Structures.  The  frequency  chosen  on  the  basis  of 
convenience  is  1.5  GHz.  Their  dimensions  are  shown  in  the 
figure  below: 
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Figure  66 
Dimensions  of  the  Wire  Crosses 

D.  EXPERIMENTAL  OBJECTIVES 

The  experiments  were  to  proceed  through  the  three  basic 
forms  (dielectric,  metal  and  wire  structures)  to  determine 

(1)  Whether  an  image  is  visible 

(2)  What  degree  of  difficulty  exists  in  obtaining  it. 

E.  EXPERIMENTAL  PROCEDURES 

The  success  or  failure  in  determining  the  first 
objective  was  seen  to  depend  principally  on  the  amount  of 
energy  transferred  to  the  element  by  the  field  and  on  its 
ability  to  reradiate  that  energy.  The  ability  of  the 
structure  to  accept  the  field  energy  rests  principally  in 
its  resonance  characteristics,  and  next  on  its  ohmic 
conductivity.  The  ability  of  the  structure  to  reradiate  the 
energy  then  depends  on  its  thermal  conductivity  and 
emissivity.  Since  ohmic  and  thermal  conductivity 
essentially  go  hand-in-hand  it  is  seen  that  there  exists 
essentially  a  three  dimensional  matrix  relationship  among 
the  parameters,  as  seen  in  Figure  67. 
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Figure  67 
Matrix  of  the  Relationship  Among  the  Parameters 

Therefore,  the  experiments  are  seen  to  be  conducted  in  a 
matrix  fashion,  choosing  an  element  with  successively  more 
complex  (less  simple)  resonance,  and  for  each  plane  level  of 
resonance  complexity  investigating  the  response  with  various 
coordinates  of  conductivity  and  emissivity. 

To  implement  this  plan  it  was  necessary  to  modify  the 
emissivity  of  the  metal  structures.  The  modifications  were 
irreversible  in  the  short  time  available,  so  a  definite 
procedure  was  followed  as  outlined  here. 

1 .   1  ._5  X M on 0£oles 

These  provide  the  simplest  and  most  predictable 
resonance  characteristics.  These  were  irradiated  at  939.39 
MHz  and  investigated  for  the  surface  thermal  patterns  in 
this  order: 

a.  The  dielectric  base  monopole. 

b.  The  brass  monopole  without  modification 

c.  The  brass  monopole  sprayed  with  flat  black 
paint  on  one  half. 
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d.  The  brass  monopole  painted  with  carbon  paint  on 
the  other  half. 

e.  The  brass  monopole  with  large  carbon  paper 
sheet  behind  it. 

2-   Aircraft  Models 

These  structures  have  significantly  greater 
resonance  complexity,  but  are  fundamental  to  the  research 
for  which  the  use  of  the  AGA  device  is  being  considered.  The 
models  used  can  be  approximated  by  simpler  structures  for 
which  there  exist  known  measurements,  notably  the  Case  1A 
and  Case  2  models  discussed  above.  These  were  investigated 
in  this  order: 

a.  Plastic  Model,  nose  down  over  ground  plane. 

b.  Plastic  Model,  tail  down  over  ground  plane. 

c.  Copper  Model  with  oxidized  surface,  nose  down, 
touching  the  ground  plane. 

d.  Copper  Model  with  oxidized  surface  in  flight 
configuration   at   varying   heights   above   the  ground  plane 

(  */?>#  Vs.  /  ^"X  /  and  other  intermediate  positions). 

e.  Copper  .Model  in  same  configurations  as 
E.  2.  b  and  c  above  but  painted  to  obtain  higher  emissivity. 
Choice  of  emissive  paint  to  be  based  on  results  of  paragraph 

E.  1  above. 

e.   Painted  Copper   Model  in  configuration  E.  2  b, 
nose  down  touching  the  ground  plane,  with  carbon  paper  sheet 
immediately  next  to  it. 
3 •   Wire  Crosses 

These  were  investigated  with  the  carbon  paper  in 
close  proximity. 

F.  EXPERIMENTAL  RESULTS 
1 .   General 

The  experiment  was  successful  in  showing  the  ease 
with  which  the  AGA  device  might  be  used  to  observe  the 
fields  of  irregular  structure.  Further  it  showed  the  carbon 
paint  and  carbon  paper  to  be  effective  transducer  media  in 
that  it   is   evident   that   they   transduce   fields   without 
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distorting  them  from  what  has  been  previously  measured  to  be 
their  correct  form.  It  is  considered  important  that  the 
relative  simplicity  of  the  method  be  appreciated.  There  is 
no  need  for  prechilling  of  the  transducer  medium.  There  is 
no  need  to  get  extremely  close  to  the  radiator,  thereby 
distorting  the  fields  or  incurring  the  effects  of  nonplanar 
wavefronts.  The  experimenter  is  not  constrained  in  source 
frequency  selection  since  there  are  no  size  limitations.  It 
is  also  evident  that  the  power  limitations  cited  for  the  DTM 
process  are  not  an  issue.    The   principal   issues   are   how 

(1)  to   optimize   the   selection   of  a  transducer  medium  and 

(2)  optimize  the  selection  of  color  assignment  to  the  sensor 
temperature  ranges  discussed  below.  Of  these  two  issues  the 
latter  is  eminently  elegible  for  subjective  manipulation. 
The  former  may  be  approached  on  a  methodical  scientific 
basis. 

2.   Guide^to, Interpretation  of  Results 

The  photographs  which  accompany  the  following 
discussion  are  conventional  Polaroid  color  photographs  of 
the  screen  of  the  AGA  device  display  monitor.  The  color 
bars  across  the  bottoms  of  the  photographs  show  the  colors 
associated  with  a  particular  value  of  the  temperature.  The 
temperature  range  is  in  discrete  .2° steps  centered  on  a 
continuously  variable  temperature  level.  The  choica  of 
color  for  a  particular  temperature  range  was  based  simply  on 
what  combination  of  colors  gave  a  percieved  clarity  of 
image.  This  assignment  can  be  varied  arbitarily.  What  is 
significant  is  that  the  device  sensed  the  temperature 
distribution  which  occurred  and  presented  it  in  a  manner 
which  then  could  be  manipulated  for  improved  human 
perception. 

3 •   £x p er iment a l_Para meters 

The  experiments  used  the  60°  corner  reflector 
structure  as  the  illuminating  source.  It  was  fed  with 
various  pcwer  sources  with  a  range  of  power  levels  from  5  to 
80  watts.  The  corner  reflector  was   consistently   placed   at 
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.4X  from  the  driven  element.  The  target  elements  were 
generally  a  distance  of  approximately  2X  from  the  driven 
element.  An  additional  parameter,  view  angle  (<x)  ,  emerged 
as  one  of  certain  interest.  It  is  defined  by  Figure  68  as 
the  angle  between  the  direction  of  arrival  of  the 
illuminating  wavefront  and  the  direction  along  which  the 
target  element  is  observed. 


Tdryt 


V 


Q> 


Viewable, 


Corner  r«£  lector 
g*-oor\4  plane 


AG-A  Device, 


Figure  68 
Top  View  of  Experimental  Setup  (Not  to  Scale) 

4 .  .  1.  5X  Monopoles 

a.   Dielectric  base  monopole 

Figure  69  shows  the  image  obtained  of  the 
dielectric  base,  relatively  fat,  1.5X  monopole.  The  view 
angle  for  this  image  is  approximately  45°.  It  shows  the 
current   distribution   of   Figure  62.   Additionally  it  shows 
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the  spillover  which  occurs  at  the  top  of  the  monopole  when 
the  current,  no  longer  able  to  continue  longitudinally  along 
the  vertical  axis  becomes  transverse  and  flows  around  the 
side  of  the  monopole  away  form  the  incident  wavefront.  This 
phenomenon  matches  that  measured  by  Burton  and  King  for  this 
same  identical  case. 


Figure  69 
AGA  Thermovision  Image  of  a  1.5X  Fat  Monopole 

b.   Brass  Monopole 

No  image  was  obtained  of  the  currents  on  the 
brass  monopole.  This  lack  of  success  is  not  considered 
conclusive  however.  Resolvable  difficulties  discussed  below 
were  incurred  with  the  carbon  paint  coating  and  with 
obtaining  a  carbon  paper  medium  large  enough  for  the  device. 
The  time  constraints  did  not  allow  these  issues  to  be 
immediately  resolved  in  this  case. 
5 .   Aircraft  Models 

a.   Plastic  Aircraft  Model 

Figures  70  and  71  show  the  Image  obtained  of 
the  plastic  aircraft  model.  Figure  70  shows  a  view  angle  of 
90°.   The   stronger   currents   apparent   on   the   port   wing 
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indicates  its  closer  proximity  to  the  illuminator.  The 
overall  pattern,  with  an  apparent  current  maximum  at  the 
Case  1A  junction,  located  in  line  with  the  trailing  edge  of 
the  wings,  is  consistent  with  the  expectations  of 
C.  3.  b.  (1)  above.  A  Case  2  image  was  not  attempted  due  to 
time  constraints. 


Figure  70 
Plastic  Aircraft  Model,  Case  1A  View  Angle  =  90° 


Figure  71 
Plastic  Aircraft  Model,  Case  1B  View  Angle 


=  90° 
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Figure  71  shows  the  aircraft  illuminated  from  its  underside. 
The  hot  spot  on  the  fuselage  next  to  the  cockpit  indicates 
the  lowest  current  maximum  of  the  Case  1A  configuration. 
Additionally  the  junction  current  maximum  is  visible  at  the 
trailing  edge  of  the  wings.  Of  possible  additional  interest 
are  the  patterns  on  the  vertical  stabilizer  indicating 
current  distributions  thereon. 

fc.   Copper  Aircraft  Model 

Figure  72  and  73  show  the  images  obtained  of 
the  copper  aircraft  model.  The  highly  oxidized  surface  of 
the  model  apparently  provided  the  high  emissivity  observed, 
while  the  thinner  metal  thickness  inhibited  flow-away  of  the 
thermal  fields.  These  results  were  particularly  encouraging 
since  the  power  levels  needed  were  much  lower  than  those 
required  for  the  plastic  model. 

Of  particular  interest  are  the  patterns 
observed  on  the  wing  surfaces  where  currents  are  not 
constrained  by  narrow  structures.  Apparent  hot  spots  appear 
near  the  trailing  edge  of  the  port  wing  as  seen  in  Figure 
72. 


Figure  72 
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Figure  73 

6 •   Wire  crosses 

Figures  74  and  75  show  the  Case  1A  and  1B  wire 
crosses.  In  all  previous  images  the  currents  on  the 
structure  were  being  observed  by  their  being  transduced  into 
a  thermal  field  on  the  structure.  Here  the  currents  are  in 
a  piece  of  carbon  paper  placed  next  to  the  cross.  The 
currents  are  apparently  induced  parallel  to  the  lines  of 
force  (LOF)  of  the  electric  field  so  that  the  image  obtained 
represents  the  charges  on  the  structure.  Figure  74  is 
completely  consistent  with  Figure  59  and  [2].  Figure  75  is 
completely  consistent  with  Figure  60  and  [2].  It  is 
significant  to  note  the  diminished  magnitude  of  the  charge 
at  the  top  of  the  Case  1A  cross,  caused  by  the  sharing  of 
charge  with  the  cross  arm  tips.  This  phenomenon  is  cited  in 
[2].  It  is  interesting  to  note  the  shift  of  charge  away 
from  the  junction  seen  in  Figure  75.  This  too  is  consistent 
with  the  measurements  of  [2]. 
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Figure  74 
Case  1A  Wire  Cross  with  Carbon  Paper  View  Angle  =  0° 


PilP 


Figure  75 
Case  1B  Wire  Cross  with  Carbon  Paper  View  Angle 


=  0° 
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8-   Techno  log  ical__  Com  plications 

Two   items  presented  relatively  minor  problems. 

a.  Emissivity  Conductive  Coating 

Yet  to  be  resolved  is  the  problem  of  an  optimum 
coating  for  either  dielectric  or  metal  base  objects.  The 
material  used  contained  acetone  as  a  base,  which  eroded  the 
polystyrene  of  the  aircraft  model.  When  placed  on  glass  or 
metal  this  paint  cracked  while  drying.  One  possible 
solution,  that  used  on  the  dielectric  base  monopole,  is  to 
first  coat  the  object  with  a  standard  lacquer66  allowing  it 
to  dry  before  being  painted  with  Eccocoat.  Aguadag  was  also 
suggested,  but  not  tried,  as  an  alternative. 

b.  Large  carbon  paper 

When  using  carbon  paper  to  sense  the  fields  in 
the  vicinity  of  the  structure  one  must  allow  indirection  of 
currents  in  the  paper  along  the  lines  of  force  of  the 
electric  field.  To  avoid  distortion  of  such  currents  the 
carbon  paper  must  be  large  enough  to  provide  a  complete 
circuit  for  the  currents  involved.  The  sheets  used  with  the 
Case  1A  cross  were  approximately  twice  as  large  as  the 
accompanying  structure.  The  sheet  used  with  the  Case  1B  was 
only  slightly  larger  than  the  cross,  while  that  used  with 
the  brass  monopole  was  approximately  the  same  size.  It 
appears  that  there  is  a  correlation  between  relative  size 
and  success  in  imaging.  It  is  suggested  that  that  success 
is  related  to  the  freedom  with  which  currents  can  be  induced 
in  the  carbon  paper. 

G.    CONCLUSION 

The  experiments  described  here  clearly  demonstrate  the 
feasibility  and  simplicity  of  using  the  AGA  Thermovision 
device  for  imaging  electromagnetic  fields.  Its  advantages 
over  DTK  are: 


66The  laquer  used  here  was  Tuf-On  747-S  Moisture  and 
Fundus  Resistant  Varnish  manufactured  by  Brooklyn  Paint  and 
Varnish  Co.,  Inc.,   50  Jay  Street,  BrooKlynP  N.  Y. 
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1 .  Clarity  _of  image 

Because  the  colors  may  be  manipulated  at  will,  the 
image  may  be  subjectively  manipulated  to  obtain  optimal 
perceived  clarity. 

2.  Ease  of  calibration 

Images  may  be  easily  measured  to  obtain  values  of 
field  intensity.  The  accuracy  of  these  measurements  is 
limited  by  step  size  of  the  temperature  increment  and  the 
calibration    of   the    device,    but    is    identifiable. 

3 .  Ability to Erasure charges or current o  n^ the 

structure 

The  methods  described  here  permit  observation  of 
current  on  the  structure  directly  or  of  charges  by  inference 
from  the  near  electric  field  observed  by  the  currents  it 
induces  in  the  carbon  paper.  DIM  is  limited  to  observing 
the  charges  on  the  structure. 
*•  •   Ease  of  observation 

Once  the  issues  of  F. 1  above  are  resolved,  the 
investigator  simply  observes- the  pattern.  There  is  no  need 
for  painstaking  presensitization  or  cooling.  Since  there 
are  no  freguency  limitations,  he  can  select  the  frequency 
most  desirable  for  his  objective.  Since  there  is  an 
apparently  lower  power  threshold  he  can  operate  at  safer 
power  levels  and  at  distances  from  the  illuminator  which 
minimize  interactions  with  the  source. 

In  conclusion,  the  use  of  the  AGA  Thermovision 
device  in  the  manner  described  here  shows  significant 
potential  as  a  means  of  observing  the  charges  and  currents 
on  structures  and  the  resulting  electromagnetic  fields. 
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